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SYBOKIS 

The  formulas  for  membrane  forces  and  displacement  in  the  parabolical 
dome  roof  and  the  conical  shell  wall  Which  are  loaded  symmetrically  with 
respect  to  the  axis  are  derived,  and  the  differential  equations  of  bend- 
ing for  the  parabolical  dome  roof  are  simplified  from  the  general  equa- 
tion of  shells  which  was  derived  by  Timosenko  .  While  the  equation  for 

the  conical  shell  wall  is  derived  from  the  concept  of  the  beams  on 

3 
elastic  foundation  . 

These  equations  arc  difficult  to  solve,  hence  the  finite  difference 
procedure  is  used*.   Thus  the  problem  is  reduced  to  the  simple  task  of 
solving  a  system  of  simultaneous  linear  algebriac  equations.  The  numeri- 
cal computation  involved  in  the  procedure  is  conside-rably  simplified  by 
two  devices.  First,  the  number  of  equations  necessary  to  attain  suffi- 
cient accuracy  is  reduced  by  an  evaluation  of  the  error  introduced  in 
substituting  central  differences  for  derivatives.  Secondly,  the  solution 
of  simultaneous  equations  is  determinated  by  using  the  digital  computer. 

The  consistent  correction  forces  at  the  edge  of  the  dome,  ring,  and 
the  ccnical  shell  wall  are  computed  according  to  the  compatibility  equa- 
tions. By  superposition  of  the  forces  found  by  the  membrane  and  bending 
theory,  the  total  forces  acting  on  the  shell  can  be  obtained. 


-'■  Numbers  on  the  upper  right  corner  of  the  sentences  refer  to  refer- 
ence listed  in  bibliography. 


INTRODUCTION 

During  the  last  twenty  years  the  shell  structure  has  achieved  extra- 
ordinary practical  importance.   The  main  reason  is  not  for  beautiful  forms 
but  for  the  chai-acteristic  interplay  of  force  in  spatial  surface  sturctures, 
which  results  in  a  considerable  saving  in  building  cost. 

The  done  is  one  kind  of  shell.   Many  massive  domes,  from  those  of  the 
pantheon  (Fig.  1)  and  St's  Peter's  (Fig.  2)  to  those  of  the  auditorium  in 
the  university  of  Illinois  (Fig.  3)  were  built  as  shells.   If  the  menthod 
of  raising  a  done  with  a  balloon  as  the  form  work  is  successfully  developed, 
the  choice  of  a  dome  roof  and  floor  will  be  the  most  economical  structure 
in  building  construction. 

Korden  dome  thickness  is  small  compared  with  the  other  dimensions. 
The  laws  governing  this  interplay  of  forces  connot  be  explained  by  the  ele- 
mentary single  dimension  stresses  analysis  of  linear  members,  mathematically 
elaborate  shell  theory  has  been  developed.   However  because  of  the  some 
mathematical  difficulties,  the  practical  design  of  these  shells  has  been 
only  based  on  the  assumption  of  a  membrane  state  of  stress  in  the  shell. 

Timoshenko  developed  the  general  equations  for  an  axisymmetric  shell 
which  carries  no  surface  loading,  but  only  edge  moments  and  shears.   These 
equations  are  homogeneous.   To  solve  a  shell  problem  a  membrane  solution 
has  to  be  superimposed  upon  the  solution  to  the  homogeneous  equations  and 
the  constant  of  integration  are  adjusted  to  suit  the  boundary  condition. 
But  the  general  equation  of  the  shell  is  difficulty  to  integrate.  Approxi- 
mate simplified  analysis  only  is  derived  for  the  spherical  dome  and  circular 
cylinder  of  constant  thickness.   Therefore  the  finite  difference  will  be 
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Fig".    2      S.    Peter,    Rome,   AD   1506-1626 


Fig.    3      The  auditorium  of   the    Uni.    of    111. 


applied  to  solve  this  particular  form  of  the  donie  roof  and  shell  wall  in 
this  report. 

Tne  dome  structure  which  is  proposed  in  this  report  is  a  parabolical 
dome  roof  with  a  conical  shell  rail  (Fig.  4).  These  two  shells  joined  by 
a  tension  ring.  The  top  radical  partterned  dome,  lite  a  cover  on  a  soup 
tureen.   This  type  of  form  can  be  used  as  gymnasium,  Auditorium,  theatre, 
exhibition  building,  etc.   The  slop  of  the  wall  can  be  constructed  as  step 
down  seating. 

Both  part  of  the  shells  are  constant  in  thickness;  it  is  six  inches 
for  the  dome  roof  and  eight  inches  for  the  shell  wall,  with  two  hundred  feet 
in  free-span  diameter.  The  ring  encircling  the  edge  of  the  dome  and  wall  is 
provided  to  take  the  horizontal  thrust;  the  conical  shell  wall  is  assumed  to 
be  fixed  on  the  top  of  the  foundation.  For  ease  in  checking,  numerical  cal- 
culations follow  the  formula  derived  in  each  part. 


DEFINITIONS  AND  GENERAL  MEMBRANE  THEORY 
FOR  SHELLS  OF  REVOLUTION 


Definitions 

A  thin  shell  is  a  curved  slab,  its  thickness  is  small  in  comparison 
with  the  other  dimensions  of  the  shell  and  with  its  radii  of  curvature,  the 
surface  thac  bisects  the  thickness  of  plate  Is  called  the  middle  surface. 

Domes  are  defined  as  thin  shells  in  the  form  of  surfaces  of  revolution. 
The  parabolical  dor.ie  roof  and  the  conical  wall  in  this  report,  the  surface 
is  described  by  revolving  an  arc  of  a  circle.   The  center  of  the  circle  is 
on  the  axis  of  rotation. 

Membrane  theory  for  shells  of  revolution 

Basic  assumption 

1.  Bending  of  the  shell  is  negligible. 

2,  Middle  surface  of  the  shell  can  be  assumed  to  suffer  only  extension, 
and  a  pure  membrane  state  of  stress  exist.   (Shearing  stresses  can  be  neglected) 

3,  Points  on  a  normal  to  the  middle  surface  before  the  deformation  sh^ll 
be  on  a  straight  line  after  the  deformation  has  taken  place  and  be  normal  to 
the  deformed  middle  surface. 

4.  Deformations  are  small  compared  to  the  shell  thickness. 

Consider  a  shell  of  small  thickness  t,  in  the  form  of  a  surface  of  re- 
volution about  the  vertical  axis.   Consider  the  equilibrium  of  a  small  element 
ds x  ds.  Suppose  P  is  the  intensity  of  loading  normal  to  the  surface  of  the 
element  of  membrane  in  the  direction  as  shown  in  Fig.  1-1. 

Then  equating  the  force  acting  on  an  element  ds x  ds  perpendicular  to  the 
surface  (Z-axis),  gives 


N'sin0/2ds    +  N*  Gin  0/2  Sin  Si  ds     +  P    ds  ds    =0  (1.1) 

0  6  z 


For     ds   =  i-    d0    =  r.dG       and     Sind0/2   =  d0  ,      Slnd6/2  =  d9       when  d6/2,   d0/2 

are  very  small. 

Hence   equation   (1.1)   becomes 


t    ,     5 


p       =     0 


(1.2) 


Meridian  force  I 


Hoop  force  N. 


Ca) 


0    0 


Fig.  1-1  Shell  of  revolution  (dome) 

Now  consider  the  vertical  equilibrium  at  section  A-A.   If  R  Is  the  total 
downward  loading  on  the  shell  above  A-A,  then 


R  =  N".2rrr„Sin0 
p    0 


-  R 


0   '  2nrQSinf 


(1.3) 


The  two  equations  can  be  solved  for  the  membrane  forces  N,  and  N  . 

0      0 

Sign  convention  adopted  is  as  follows: 

For  N'  and  N*  Tension  + 

<t  a 

Compression 

For  r_,  r         Concave  inward   f 

Concave  downward 

For  P  Acting  outward  + 

z 

Acting  inward 

Displacement  in  syiantrically  loaded  shells  having  the 
form  of  a  surface  of  revolution 

In  the  case  of  symmetrical  deformation  of  a  shell,  a  small  displacement 
of  a  point  can  be  resolved  into  two  components  v  in  the  direction  of  the 
tangent  to  the  meridian,  and  w  in  the  direction  of  the  normal  to  the  middle 
surface  (Fig.  1-2).   The  change  in  the  clement  due  to  the  difference  in  radi- 
cal displacements  of  the  points  A  and  B  can  be  neglected  as  a  small  quantity 
of  higher  order.   Thus  the  change  in  length  of  the  element  AB  due  to  deforma- 
tion is 

—  dp"  -  w  it 
d0 

Therefore   the   strain  of  the   shell    in  the  meridional  direction   is: 


P  dp  1  r  c\0       r. 


The  radius  r  of  the  circle  increases  by  the  amount 


v  Cos  $   -  wSin  I 


Hence 
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£        =     — (vCos  0  -  wSin0) 

0         ro 


(1.5) 


Or  substituting     r.    =  r„Sln0 


C .      =     —  Cot0    - 
8  r„ 


(1.6) 


Eliminating  w  from  Eqs.    (1.4)   &   (1.6) 


dv 


vCot0     =     rxt(-   r2   G6 


0.7) 


Fig.    1-2     D is placement  of   the   shell 

The  strain  components   6.  and  6     can  be   expressed    in  terms   of   the   forces   N, 
0  9  0 

and   N*    by  applying   Hook's   law.      This   gives 


s  =  £(n0-'n;> 


£<■;  ■%■> 


(1.8) 
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Where  i    is   the    polsson's   ratio.      Substituting    in  Eq.     (1.7) 


civ 

-—  -  v  Cot  </.>     -- 

tl0 


i[^(r,    +*r2)-«J(r2+Vr1>]  (1.9) 


Denoting  the  right-hand  side  of  this  equation  by  £(0),  hence  yields 


~   -  vCol:0  =  £(0) 
d0 


The   general  solution  will   be 


«.#[;£§.«  +  .] 


(1.10) 


whore  c  is  a  constant  to  bo  determined  by  the  support  conditions.   From 
Eq.  (1.6) 


w  =  vCot  0  -   r  6 


'Cot  (8  -  —  (K^  -t/Nj)  (1.11) 


The  meridian  rotation  A0  can  be  expressed  in  terms  of  displacement 


A0  =  ■*-  +  -£, 
ri   V" 


(1.12) 


The  horizonal  raovemsnt  <&H    can   be  derived  directly  from  Eq.    (1.8) 

(1.13) 
The  meridian  rotation  at   the  edge   will    be,    from  Eq.    (1.12)   with     v  -  0 


r    Sin  {i 


A0    « 


dw  Cot  0  dv  d      r    9 


r.dp 


r,       d/,         r 


(1.14) 
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From  Eq.    (1.9),  with     v  =  0 


%    =    &>;<*!  +Vr2>-»J(r2Wrl>]  (!.»« 


Substitute  Eq.  (1.15)  into  Eq.  (1.14),  hence  gives 


A*  -  St^i  ♦V*  -  «;(r2  +Vri)3  -  --^C^cn;  -y#]    a.») 


Where  only  the  horizontal  movement  is  required,  it  is  only  necessary  to 
compute  N.  and  N  at  the  edge. 

With  these  Eqs,  (1.10),  (1.11),  (1.12),  (1.13),  (1.16)  the  displacement 
of  shells  due  to  membrane  theory  can  be  solved. 
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MEMBRANE  AND  RIGOROUS  ANALYSIS  FOR  THE  PARABOLICAL  DOME 

« 
The  property  o f  the  dome 

The  upper  part  of  the  dome  is  parabolical  shape  which  is  shovm  in  the 
following  figure 


Fig,  ?.-!   Sect,  through  the  dome 


Assume  the  equation  y  =  kx   Is  the  function  of  the  curve  of  the  dome,  r   is 
the  radius  of  curvature,  r   is  the  radius  of  curvature  crossed  the  axis  of 
revolution  of  the  dome. 


dy 


=  2  k  x  =  Tan  t      =2  /Ty 


dx2 


=  2  k 
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S!n(5 


,;>■■ 


dy 


dy/dx 


V  dy2  +  dx2     V  (dy/dx)  +  1 


Tan(S 


2V  ky 


V  Tan20  +  1     /**?  +  X 


(2.1) 


Cos0   = 


7  1  4  Tan' 


2.     V  4ky  +  1 


(2.2; 


c*  ♦€»"] 


3/2 


3/2        3 
(1  +  4ky)        Sec  0 


a7 

.  2 
dx 


2k 


2k 


(2.3) 


,dx 


Sec  t 


'2    =  >V  *  H  <$     *   J  i ,/  »  +  TRy  -   ^v1  «  4^  =  —   ».*) 


Assume  the  thickness  of  the  dome  Is  6  Inches,  pnd  the  radius  of  revolution 
at  base  is  100' 

Let  the  curve  of  the  dome  to  be  divided  into  10  divisions,  the  property  of 
each  point  is  calculated  as  shown  in  table  2.1. 


4   3   2   1 


Fig.  2-2.  The  division  of  the  dome  in  forces  analysis 
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TABLE  2.1  PRCPCRTY  CF  THE  DOME 

RAD        SINX         COSX  SECX  TANX       DEG 

.OOOOOOOO   0.00000000   1.00<  0   0  1.00000000   O.COOOOOOO   0.00 

.'.5235988    .05233596    .99862950  1.00137240  .05240778   3.00 

.1<.  471976    .10452847    .99452167  l.uu550630  .10510424   6.00 

.15'/G7963    .15643446    .98768631  1.01246520  .15838444   9.00 

.20943950    .20791168    .97814760  1.U2234C60  .21255655  12.00 

.26179937    .25881903    .96592585  1.03527620  .26794917  15.00 

.314]t-924    .30901697    .95105652  1.1-5146220  .32491-367  18.00 

.36651911    .35836792    .93358  42  1.07114500  .38386401  21.00 

.41687698    .40673661    .91354548  1.O9463630  .44522664  24.00 

.47123685    .'5399046    .89100654  1.12232620  .50952540  27.00 

.52359672    .49999996    .86602542  1.15470050  .57735021  30.00 

DEG        P]               R2  X  Y 

.'  '   173.20510000   173.20510  0.  0. OOOOOOOO 

3.o0   173.91921o00   173. 44281000  9.07729530  .23786046 

6.wU   176.o83u9l.00   174.15917000  lo.2045920u  .95668998 

9.0(   179.76329000   175.36414000  27.43299500  2.17247990 

12. rr       185.0748900C   177.07461000  36.81588000  3.91272830 

15.00   192.18937000   179.31512  00  46.41016500  6.21778280 

18.00   201.34539000   182.11862000  56.2777440U  9.14287300 

21.00   212.8656o>.00   165.52778000  66.467205O0  12.76102300 

24.00   227.18006000   189.59659000  77.115«740u  17.16709800 

27.00   244.86007000   194.39262  00  80.25239500  22.48341700 

30.00       266.6666700C        200.    00    2    00  100.1        01      00  28.86751000 


16 


Membrane  forces  due  to  dead  and  live  load 


let  the  uniform  load  acting  on  the  dome  as  shown  in  Fig.  2-3 


Fig.  2-3   Uniform  load  acting  on  the  dome 


Loading  q 

Loading  in  meridional  direction  P.  =  qSin0 

End  load  at  each  section  is  R.   S  is  the  length  of  the  curve. 

For 


y  =  kx  , 


dy  =  2  k  x  dx  , 


2     2  2  2 

ds  =  J  dx  +  dy   =  7  1  +  4k  x  dx 
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Hence 


,    2   2 
1-iAk  x 


2   2 

R     =     /  2rr  q  x  v  1   +  4k  x     dx 

0 


^i-  [(1   *  4kV)3/2   -   1]  (2.5) 

6  k 


The  membrane   force    in  meridian  direction   (Eq.    (1.3)) 


0  2TTroSin0  4nkx 


The  membrane   force    In  the  direction  tangent   to  the  circular  cross   section 
can  be  derived  from  Eq.    (1.2) 

■   R  Vo 


Zrrr   Sin  $      sil& 


2nr    Sln2(J 


q  X  Cot0 


4y  /l+  4ky       2  k 


(2.7) 


For  dead   load   75    lbs/sq.    ft.    and   live   load   30   lbs/sq.ft,   use   these  Eqs.    (2.6), 
(2.5)   and    (2.7)    to  calculate   the  stresses  resultants  which   is   given    in  table 
2.2   and   2.3. 
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TABLE  2.2  MEMBRANE  FCRCEb  DUE  TC  D.L. 


DEG 

R 

M* 

'Jfl 

.  1 

O.OOOi 

-6495. 1913 

-6495.1913 

3.10 

19427.2990 

-6508.4244 

-6499.7860 

6.0C 

78301.0750 

-6548.9593 

-6512.9790 

9  .  C  ( 

178424. 2600 

-6617.1088 

-6535.2060 

12.  0C 

322936.8000 

-6714.6510 

-6565.9860 

1  5  .  C  0 

516500.5800 

-6643.5508 

-6605.2620 

18.00 

765609. 3000 

-7i;06.6147 

-6652.8390 

21.00 

1079041.2000 

-7207.6040 

-6708.4340 

24.00 

1468496.2000 

-7451. .3656 

-6771.7320 

27.01' 

1949518.1000 

-7744.1581 

-6842.3530 

3  >.<  i 

2542808.0000 

-8094.0098 

-6919.8750 

TABLE  2.3  MEMBRANE  FCRCEi  DUE  TC  L.L. 


DEG 

R 

Nf 

N# 

. 

O.OOOC 

-2598. u765 

-2598.076? 

3.0( 

7770.9199 

-2603.3698 

-2599.9141 

6.00 

31320.4300 

-2619.5637 

-2605.1915 

9.0( 

71369. 715C 

-2646.8436 

-2614.0821 

12.0  0 

129174.72  00 

-2685.8604 

-2626.3940 

15.00 

2C  66.  0.2300 

-2737. 42u3 

-2642.1048 

3  8.0  0 

306243.7300 

-2802.6459 

-2661.1351 

:   1  .   :. 

431616. 4a00 

-2883. u415 

-2683.3736 

2  4  .  0  C 

587398.5000 

-2980. 5465 

-2708.6924 

'7.  ■ 

77981  7.2700 

-3097.6634 

-2736.9409 

3  .00 

I1  17123.2000 

-3237.6040 

-2767.9499 

Membrane  forces  due  to  snow  load 


The  snow  load  P  acting  on  the  dome  Is  shown  In  Fig.  2-4. 
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For 


P   =  P  Cos  0  , 
z 

R  =  PttX 


P   =  PSin  0  Cos  0  , 


(2.8) 


From  Eqs.  (1.3)  and  (1.2),  hence 


Pttx 


Px 


0      2rrroSin0      2nr,,Sin2 


2r  Sin"0 


(2.9) 


2TTriSin20    Z  Sinf     2Trr1Sin20    2 


r.PCos  0  (2.10) 


For      P  =  30   Ibs/sq.ft,    the   solution  of   the  above  Eqs.    is   shown  in  table  2.4. 


4> 


Qnnmrc 


■Pis^ 


P=30  lbs/ft 

ZLXXxxrxxo 


Fig.    2-4     Snow  load   acting  on  the  doi-e 
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TABLE  2.h    UEN"DRANE  FORCES  DUE  TC  SNGW  LCAD 


DEG 

R 

♦ 

N« 

. 

.  ■ 

-2598.0765 

-2598.0765 

7765.7617 

-2601.6419 

-2594.5165 

■  . 

31234. 

-2612. Jc  ^6 

-2583.8437 

9.00 

70927.9780 

-2630.4622 

-2565.0902 

1?.  ( 

127744.2900 

-2656.1191 

-2541.3022 

203UL0.6200 

-266v. 7268 

-2509.5495 

It. 

29851.0.1100 

-2731. 7794 

-2470.9174 

21. CO 

416626. £ 

-2 762.9168 

-2425.5135 

24.1 

560478.1700 

-2843.9489 

-2373.4613 

2  7.00 

734(  /,7.4C  0 

-2916.0893 

-2314.9030 

30.  rr 

942477.8100 

-3000. 

-2250.0002 
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Displacement  from  the  membrane  theory 
From  Eq.  (1.13)  the  horizontal  displacement 
r  Sin0 

AH  -  -ir- (»;  -Vty  (2.11) 

And  from  Eq.  (1.16),  the  angle  rotation  in  meridian  direction 


For 


-4*(~!4-)  = 3— ~  [(Sec0Csc20  ..  Cso20Cos2Cl  -  3  Sec  0 

rid0  Sln0      ri  l2k2Et  d0 

-^(See  0Csc0  -  Sec  0  Csc  0  Cot  0)] 


[  (-  2  Csc  0  +  Csc  0  Sec  0  +  2  Cos  0  Csc  0  +  2  Cot  0 


12r  k2Et 


-  3  Sec  0  Tan  0)  -  V  (-  Sec  0  Csc  0  Cot  0  +  3  Sec  *0 


3  2 

+  Csc  0  Sec  0  -  Csc  0  Sec  0  +  Csc  0  Cot  0)] 


Let 


Hence 


r  d0vsin0 


A0   =  ?^r^(ri  +^r2)  "  Ne(r2  +Vri^  • A  (2-13) 

The  sign  convention  soe  Fig.  2-5.   For  the  numerical  calculation  of  Eq.  (2.11) 
and  (2.13)  see  table  2.5  and  2.6. 
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on:ied    position  of 
membrane  condition 
under  gravity   loc.d 

(a) 


Cb) 


(c) 


Fig.  2-5  Sign  convention  of  the  deformation  correspond 
to  the  edge  forces 
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T4BLE  2.5  DISPLACEMENT  H«CM  MEMBRANE  THEORY  (u.L.) 


DFG  E-AH 
0  .  C  ( 

3.  -94369.405000 

6.'   li  -189443.800000 

9.0C  -285949.7001  ' 

12.:  (  -384582.7900  )i 

15. 0(  -486053.470000 

|H.  -591086.970000 

21.00  -700364.68001  0 

24. 'JO  -614568. 630OOU 

27.00  -934331.880000 

3i  .01  -1060214.600^00 


0  .OooOOO 
2.615126 
-9.816866 
-15. 5717 j 5 
-2o. 425761 
-24.665^59 
-26.27t.6b7 
-31.191613 
-33.363228 
-34.621466 
-35.5^4618 


T4BLE  2.6  DISPLACEMENT  FROM  MEMBRANE  THEORY  (L.L.) 


0E&  E.Ari 

.00  O.OOOoCO 

3.01  -37747.758000 

6.  '  -75777.517000 

9.1  I  -114379.86001  0 

12.0i  -153833.09 

1  !  .  -194423.390000 

18.01  -236434. 740o00 

21.00  -280145.870000 

24.00  -325827.3901/00 
27. uO  -373732.710000 

30.1  1  -424085.820000 


A«f 

o.O^  0O0Q 

5.854326 

-3.331127 

-6.053862 

-8.096722 

-9. ^27795 

-11.285249 

-12.460994 

-13.341567 

-13.919031 

-14.193610 
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Differential  equations  of  banding  in  the  dome 

For  the  loads  symmetrical  about  the  axis  of  rotation,  the  formula  was 
derived  by  Timoshenko 


r2  d2U    1  r  d  V  ,  r2  _  ,  ,  r2  dt  dU 


J-[!LCot20.V-4l!lcot0]u  f  ^2f 
fj L  r2  h  d0     J      fj  <UI 

dr-  Vr2  dt    , 


+  -—[V  — i  +Vr.  Cot  0  +  r.  Cot  0 1  K. 

r1LV  d0      2    "    1    "    t  d0J  0 

r,  dN'    ,   dr  r  dt 

-  —  — 2-   -  —  [—1   +  r,  Cot0  +Vr,  Cot0  -  —  —  1  K.   =  EtV    (2.14) 
r  d0    r  Ld0    2        1       t  c!0  -'  0 


^4  +  I[i(^^cot«  +  3-!?-^ 

r2  d02    r1Ld0  rx    r^  r^t  d0 "  d0 

l_u]  _   SiDotl  dt   fi  Cofc2     =  _U 
d0   r2  D 


Angle  of  rotation  in  ineridian  direction 

Q.  is  shearing  force 
s0 

The  flexural  rigidity  for  the  shell 


Because  of  constant  thickness,  — r  -   0,  ^—  is  small  when  compared  with 

another  term,  therefore  cancel  it  in  calculation.   Also  neglect  the  effect 
of  the  membrane  stress  resultants  on  bending.   Then  all  the  term  I.".  and  N 


ri 

Where 

V 

= 

0y 

u 

= 

r2Q0 

D 

= 

Et3 

(2U  -,/2) 
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arc  dropped  from  the  equation 

2 

^"+Il[i..^+!2cot(5]^-!l[^Cot20]u      .     BtA  (2.16) 

d02       r2Lcl0rx       rx  d0       r2     r?  r? 

4  +  !i[i.(!2)+!iCot0]dv_^  ^Cot2  ]v  =  _^I 

_    .£  -I"        "-    ^(7,      f  t-  -^    Art  v        *-    -r-  -1  n    *- 


,iz        r_     dR  r_  r.  J  do1,       r.  "-  r.  J  n  r_ 

For 


d02       r2     c^  ri         ri  J  d0       r2  L  r2  D  r? 


Hence 


2             Sec  0.2k                 1  „     2. 

—     =     = =     _-     =     Cos  I 

1  Sec  0- 2k  Sec  0 


-^r  —     =     -2Cos0Sin0      =     -Sin  20 


2 

^7  +  -V  C-  Sin  2<*  +  C03^  Cot  fl  77  "  -VC-V  Cot20]  v 

d02   Cos20                      <"*  Cos20  Cos20 

=     EtVSec  0/2k  (2.18) 

2 

— *•  +  a-  [-  Sin  20  +  Cos  0  Cot  0]  — =*-[" — V"  CotVl  V 

d02        Cos^0                                                 "    df  Cos20~Cos20 


To  simplify  the  above  equation,    then  get 

2 
—  +  [-2Tan0  +  Cot0]SL    .  _JL_U      =     EtVSec50/2k  (2.20) 

d02  d0         Sin220 


2. 


5, 


— ^+[-2Tan0  +  Cot  0]  ^ 1_V      =     -£Sec_0  (2>21) 

d02      .  d0         s.n22(j  D    2k 
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The  integration  of  these  two  equations  is  difficult.   Finite  differences 
can  be  applied  to  solution  of  these  problems. 

Application  of  finite  difference 
equations  hi  bending  analysis 

The  application  of  finite  difference  equations  to  the  solution  of 
difficult  structural  problems  is  in  large  measure  comparable  to  the  technique 
now  used  to  surmount  mathematical  difficulties  in  the  solution  of  complicated 
dlf ferentii'l  equations.   Essentially,  the  technique  employed  consists  of 
replacing  the  derivatives  of  differential  equation  by  its  central  difference 
equivalent.  The  problem  is  thus  reduced  to  the  simple  task  of  solving  a 
system  of  simultaneous  linear  algebraic  equations. 

Just  as  the  replacement  of  an  integral  by  the  summation  procedure 
involves  the  use  of  average  value  of  the  ordinate,  so  the  replacement  of  a 
derivative  by  finite  differences  is  based  on  talcing  the  difference  of  average 
value  of  the  ordinate.   In  this  light  then,  it  is  evident  from  geometrical 
considerations  of  Fig,  2-6,  That  if  y  =  f(x),  then 


« 

y=f(xk 

i    v 

1 

/ 

/ 

t 

(M 

CM 
CI 

r-J 

JT 

<H 

<M 

<M 

U-\ 

%r 

4 

X-i                             ' 

i 
4 

, 

x=i+?g 

■ 

i   (yt^"W 


Fig.  6 
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Let  tho  distance  between  .each  point  in  the  x  ordinate  is  h,  then 

dx,   '     h 

In  which  y.o.  represents   the  ordinate  at     x   =  i+'j     and  '=.  means  approximately 
equal   to.      By  repeating   this   process,    it  naturally  follows   that 

il    ■      t&i      ■      yi+i  "  yi      yiJLl±k         yi-:-i  "  2yi  +  yi-i 
dx2       '    te  d*   *  h2  "       h2  h2 

dxj       '  h3 

.4 


A     .       yi42    -  %+l  +  6yi   -  %-!   4   h-1 


dx.  h 


From  these  expression  of  derivative,   the  finite  difference  equation  of 
equations    (2.20)   and    (2,21)  will   therefore  be 


°i+I   "  2Ui   +  Ui-1                                                  "■  k1    -    0,    ,  4 

-22 L. Li  +   (-2Tan0.    +  Cot0,)-iiJ bl  . -2 U. 

(A0)2                                  >                l          2*0  sta2^)     ' 

5 
=     EhV    Sec  0  /21c 


V,  ,.    -   2V.    +  V.  V  -  V 

^+(-2  1*,*,     +  Cot  CO     1+1 


(*«'  '  2    '  Sin2(2,5.) 

B     Sec  0 

D         2k 
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Collecting  similar  terms  of   U's  and  V's  and  divide  both  side  by  E 

u,  .,  1).  ■)  7 

-~[l   +  (-2Tan0.    +  Cot0j)A0/2]  -  _L[  2  +  440  /Sin    (20.)] 


«,    ,  ,         5 

+  — —  [l   +   (2Tan0j   -  Cot0.)A0/2]     -   tV.A0    Sec  0  /2k     =     0  (2.22) 


V.      [l  +  (-2  Tan0.    +  Cot0.)A0/2]    -  V.[2   +  4A02/Sin220,] 
i  -fi  ~  i  II  i 

6(1    ••  1?)&£-  Sec50. 
+  V       [l   +  2Tan0.    -   Cot0.)A0/2]   + i.  =     0  (2.23) 

l       *  Et  k 


Since  the  trigonometric  functions  for  the  various  values  of  0  can  be 
readily  evaluated,  two  difference  equations  for  each  point  can  be  written. 
Because  of  symmetry,  U  and  V  are  zero  at  0=0.   It  follows  that  the  equations 
atA0  =  0  are  superfluous,  At  boundary,  it  is  assumed  AH  =  1,  A0  =  o; 
AH  =  0,  A0  =  1  in  order  to  get  the  expression  of  stiffness  both  for  dis- 
placement and  rotation  of  the  doine  shell. 
Since 

«0  =  -  —  UCot  0  (2.24) 


_1_  d_U 
rl  d0 


r„Sin0 


(2.25) 


EAH  =  - (Hc  -V.N0)  (2.26) 


Substitute  (2.24)  and  (2.25)  into  Eq.  (2.25) 
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r„  Sin0         ,     U  -   U 

ah  .  -2 — (.JL_ta — 1=1 

Et  rx  2    0 


+  -  Cot0) 

r2 


r2Si"0  1>COs0  r2Sln0 

■U.    .    +  - -0.     +  n-r U. 


2r   AjJtE.i+1  tE         i         240  tr  E    i-1 

For    Ah  =  1,    A0  -  0,   hence 

^  Cos  0 


r„  Sin  0 


r    Sin  J 
2 

-U.     +  — - U. 


■2r,A0t    i+l  Et         1         E-2A0tr       i-1 


(2.27) 


(2.28) 


(2.29) 


Now  there  are  10  points,  and  therefore  22  simultaneous  equations  for  total 
22  unknowns  are  obtained.   Let  the  coefficient  of  Eqs,  (2.22)  and  (2.28)  to 
be  C.  and  D.  for  Eqs.  (2.23)  and  (2.29),  then  writing  in  matrix  form,  it  is: 


12  3 

cn       cu       cu 


11 


2  3 

D  D 

11  11 


'VI 

'  4  1 

C2 

V2 

4 

D2 

U3 

4 

< 

*  =  « 

\ 

Ull 

4 

Cll 

s% 

K,j 

(2.30) 


I.e. 


V<V  ■  (V 


(2.31) 
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Hence 

(U.)  =   (A   )*'(G.)  (2.32) 

Use  the  boundary  condition  AH  =  0  ,  A0  =  1  so  U's  and  V's  for  the  unit 
rotation  at  edge  of  the  dome  shell  can  be  got. 

Then  substitute  these  U's  and  V's  Into  following  equation,  the  force  due  to 
these  boundary  displacement  and  rotation  can  be  obtained 

M   r,  -  — U.  Cot0  (2.33) 

0      r2  l 

1  dU-       (U,  t,  "  U-  ,> 

(2.34) 


f,  d0         r  -2A0 


Cot0.  JdV  Cot0  > 

,     dV  V.  Cot  0  .  ,  V.Cottf 


+  V-i )   =  -D[r^r(vi+i  -  v-  ,)+V- ] 

r,,         L-2r1A0   l+l    i-l      r„ 

(2.36) 


'j  d0         j-2  "."if   "*     i-j.       i„ 


The  matrix  A,  A~  ,  U,  V  and  the  solution  of  Kqs.  (2.33)  to  (2.36)  are 
listed  in  table  2.7  to  2.13. 
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TABLE  2. T  CCEFFFCIEMT  MATRIX  OF  AX=G 

-3.0036633C  -.23905952  1. 49679890  0.00000000 

.00000000  O.OOOOOOOo  0.0(    I  o.ooocoooo 

.  :   OOl  i  U.OOOOOOOd  O.UOuOOOOO  0.00000000 

.00000001  0.00000000  O.OOuOOOOO  0.00000000 

,00000000  O.OOOOOOI  O.oOuOOOOi  0.00000000 

.00000000  0.00000000  43.03071300      -3.00366330 

1.4967989'  O.U  0.00000000 

.00000000  0. 00000000  O.OI  0.00000000 

.00000:-  O.OOOOOt  I  O.OC    000  O.OCOOOOOO 

.00000000  0.00000000  0.00000000  0.00000000 

.ooooocoo  •  .  000(  OOi  O.OOi  I  I  O.Oi  0  10000 

.7564178.  0.00000000  -2.25368760  -.24403729 

1.24358220  O.OOOOl  I  0.00000000  C. 00000000 

,00000000  0.00000000  O.OuoOoOOO  0.00000000 

.00000000  O.OOOOOOC  O.ujL'i-L'Udi,  o.oooocooo 

.00000000  0.00'jOOCOu  0.ou*jOuiJGo  0.00000000 

.0001  G.GOoOuOOO  O.uu,;0u00u  .75641760 

43.926713  I  -2.25368760  0.00000000  1.24358220 

.00000000  0.00   '000  O.Oo  OOOOI  0.00000000 

.OOOOC  0.0001  O.OOJOOOOO  0.00000000 

,0000000  0.00000000  O.OC  C  0.00000000 

.1  OOOOOOO  0.  lOOOOOOt  O.OC  '  O.OOOOOOOC 

.00000001  0. 00000000  .84299940  0.00000000 

-2.11483980  -.25259716  1.15700060  0.00000000 

,00000000  0,00000000  O.OOOOoOOO  0.00000000 


.01 

. .  . 

.  I  I  0 
.0" 
.00000  10 

.00000 
.88796250 

1 .1120375r 

.1  01 
.CO!.' 

47.7285780b 
.  I 

.000    000 i 
.  300 

.1   301  0      ' 
,(  OC 

-2.  4386490 
.000  ' 
.  DOOOO 

,00( 


.     i      OCOOt 

.   I 

0.001  ' 
.84299943 

1.1570C060 

.     .. 
. 

O.OOOl  00 

. 

0.001 

0.  oc 

.    0  100000 
-2.06628730 
.0001     o< 

.   ■  ;ooooo 

.'  0 

OCOCC 

.    ■■   000 

-.28236385 
0.     001      0 
).00    (  OOt 
3.00CI 
0.001 00000 


0  .  Oi 
0  .  01 

45.<+b748o00 
0  .  u  0  I 
0.00 

0.0  (J 
0.00 

-2. -662873L 
0.1 
0  .  0  0 1 

... 
... 

0.0 

O.OOOOOC  OC 
0. UOOOOOOU 
O.OOOl 
0.0( 
0  .01 

.9163251 
1.L836750U 
0.01  000 

0  .  '-  0 
0  .00 
0  .  0  0  ( 
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. 

"ooo 

0. 

1 

. 

,0000  .000 

2. 

,11463980 

0, 

,00000000 

, 

,00000000 

0, 

,001 

0, 

'000 

0, 

oooo 

.26515877 

0, 

,00000000 

II. 

-000 

0, 

.00000000 

' 

. 

.88796250 

1 

.11203750 

0 

.00000000 

0 

.00000000 

0 

.noooc ooo 

0 

.00000000 

0 

•ooooc 

0 

.00000000 

0 

.   i  0000 

0 

.00000000 

0 

.00000000 

0 

.00000000 

•  out-.  I  I 

.oc 

,000<  OOOi 

.9364392C 

1.06356080 

.'■  001  I 
.00000000 
.OOOOOC  10 
,0i 
5*. 92483800 
,000    . 

.<    oooooo 

.OOOl 

.        100000 
-2.02449260 

,0001 
,00000000 

.    '  I  001  00 
.000!'.''    II 

.    ■ ' OOOOG 


.91632500 
1.    8367500 

0.00       I 

rOg 

0. 

0.00000000 

o.oooooooe 

0.001  0<  ( 
O.i  3 

O.OOi      iOi 
0.00000000 
0.00    I  '  ' 
-2.03174090 
0.00000000 
...    0    0  jOO'j 
O.OOI  0    0< 
0.00000000 
0.0000000 
-•33478751 
,i  OOOOOOO 

:,.    IOI    Ol    I 

i  .  iqi  oooo 

0  .  00  OOl  Ol 

.9518980' 
1.1  4B1C201, 
.     3    Oi 
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50.82549300 

-2.04386490 

0  .  0  C 

0.00000000 

0  •  0  0  0  0  0  0  0  0 

0. 00000000 

0. 30000000 

0.00000000 

0 • oOO  0  C  0  00 

0.00000000 

0.  UUOOOOOO 

0. 00000000 

-2.03174090 

-.30513799 

O.OOX  "   . 

0.00000000 

o.oc 

".t-onnrncc 

0.00 

0.  ©0  00  OOOO 

3 .  Ol   OOOl 

0,(  OOOOOOO 

0.00000000 

.93643920 

0.  uuoOuOOo 

1.06356080 

O.OOC  ouooo 

0.00000000 

0 .  ^0  ..  QOOOQ 

0.00000000 

o.oo :  l 

0.00000000 

0.00000000 

0.00000000 

•951898C 

0.00000000 

1 .u481020u 

0.00000000 

. 

0. 00000000 

0  •  0  C 

0.00000000 

0.00  ' 

0. 00000000 

0.00 

o.oooooooc 

60.26175300 

-2.02449260 

0.  i  ■  jOOO 

0.00000000 

0.00000000 
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. I  ' OOOul 
. ,  o  o  0  0  0  i 
.OOOOf  )i 
.96451]  1 
1.  (,3548890 
.  '  C  00(  ; 
,(  00(  001  0 
,0000(  " 

67.166  25C 

.00000000 
,000( 

.  II  oooooc 

.  00 00 00 1 
-2.;  167549 
.000001  I  0 
.1  001 
.00000 
.0001  0 
.00000000 
. 0000000  0 
.001,". 

.  ( 
.001 
•  OOC  ' 


.  ( 

100000 

O.OO1. 

.   000000 
0.00  I 
0.1  01 
O.OC  i  OOi 

.  •  ■  ioo 

-2.  -1985680 
00000 

0.0000000- 
C .001  00000 
0. 00000000 
-.42278880 
0.00000000 

.    I   ' 
.  ■    oooco 

0.00000000 

.97529770 

1.02470230 

0.000(  I  i  o 
0.00(  01  (  0 
.    I      00000 


0. u0o( -000 
O.OC         i  Oi 
0.0  c 
-2.^1985680 
0  .  OOC   i< 

O.OOC 

O.Oi 

O.OC 

0.00  K    0 1  ■■ 

O.oo 

O.OOOl    )00  . 

O.oOuOOOOo 

O.uOOC 

.97529770 
1.02470230 
O.OI 
0.00    ' 
0.00    :     000 
0.1  01 

76.  Iul98400 

0.  OOi  ooooo 
0.000       I  c 
0 . o  OuO    -     - 
0.  OOOOiOOt 
0  .  0  0 1 


O.OOC 
0.0000  101 
0.00000000 
-.37314458 
0.00000000 
0. 00000000 
0.00000000 
0.00000000 
.96451110 
1.03548890 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
0.00000000 
'OOOOO 
0. 00000000 
0.00000000 
-2.01675490 
0.00000000 
'0000 
0.00000000 
0.00000000 
0.00000000 


r5 


.98486501; 
1.Q151150I 
.OO0000OC 
.00000000 
.000! 
.OOOOOOUO 

37.72981600 
.00000000 
.  00000 
.00000000 
.00000000 
7.16197200 

-7. 36197200 
,C  00  I 

.0001  i 

.(  0000000 

.OOOOOOUO 
.  00  C  i 


.  '  OPOOO 

.   ; 

o.oooooood 

0.000( 
O.OOoOOOOu 
0. OOOOOOUO 
-2.0X462160 

. 
O.i  0-000000 
0.00000000 
0.000' 

.  0 
0.0  0'. 

0.00  000000 
0.00'-'. 
0.00000000 
0. . fffOOOOO 
1.00. ' . 


-2.01462160 
O.i.i 

0  .  000  01  0  0(J 
D/.UOOOCOOC 

0. UutUuUOO 

O.uuOOUOO 0 
O.uuuOOOOO 
O.OOC     .    >( 
0.01  . 

0.1 

0.00<  ""  '.     0 
.80.        I      6 
0.0001 

. 

0. 00( 

0. UOUQOOOC 

0.000' 

0  .  o  0  u  i 


-.48738785 
0.00000000 
0.00000000 
0.00000000 
0.00000000 

.98488500 
1.C1511500 
0.00000000 
0.00000000 
O.OOOOnroo 

.<  '  OOOOO 

o.oooonnoo 

0.00000000 
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-0.45654129E- 

03 

0.53261312E- 

05 

-0.54854827E- 

05 

0.13C03242E- 

01 

0.29138411E- 

CI 

-0.59514647 

0.12736417E- 

■01 

0.34491028E- 

■02 

-0.57363860E- 

■04 

-0.64607186E- 

•04 

0.74844485E- 

-03 

-0.13468997E- 

-03 

-0.87730331E- 

-01 

0.27949766E- 

-02 

0.73911638E' 

-04 

-0.51967947E 

-05 

-0.25953006E 

-04 

0.26571037E 

-01 

0.99783302E 

-02 

-0.26789455E 

+  01 

0.33448401E 

-01 
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TABLE  2.g   E.U  AND  E.V  VECTOR.  FCRAH=1.  A+=0 

E.U  E.V 

. 00  000ul8912760  .0 0000137435850 

.00000059903112  -.00000267917660 

<)00u44480459  -.00oo2685071800 

-.0000054  8548270  -.O00O64607186  00 

-.o000259530060o  .00013682969000 

-.00  000665  774200  .00152992  0700  00 

.00044524965000  .',0324927080000 

.001905  78  860000  "..  2 0  71  338000000 

-.0041614268  -.16704788000000 

-.07892756900000  O.OOOOOOouOOOOOO 

-.1526040600  6.98327210000000 

TABLE  2.10  E.U  AND  E.V  VECTOR.  FOR  AH" 0  A<)>  =  1 

E.U  E.V 

.00000077846737  .0 000 02 55 7 79340 

.00000197068490  -.U0001724698500 

-.00  000028662721  -.00010242165000 

-.000024  32  0  53200  -.000163382  73000 

-.00008391898300  .0  0082203663000 

.00007648043200  .  .562444170000 

.00 183365 56 OOOC  .00607107300000 

.00541162840000  -.  '9880945200000 

-.  2675844100000  -.54941582000000 

-.28*50217000000  l.OCOOOOOoOOOOOO 

-.0  68537642  00000  2  7. 1053580o000000 
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TA?Lfc"    2.11    FOCES    DUE    TO    BOUNDARY    DISP.    FORAH=l 

DEG  N+  Ne 

.00  '...OO'1  O.oooCGOOO 

3.0C  -.00000002  -.00000003 

6.00  -.00000003  -.oouOOOOl 

9.0i  -.00000002  .00000032 

12.1.0  .00000015  .oOo00136 

tv.  0  .00000054  .oooooul6 

18.00  .00000015  -.00002235 

21. OC  -.00oo0625  -.00008588 

24. 0(  -.00002258  .00019364 

27.00  .000042ol  .00315242 

.00068353  .oo531571 

TABLE  2.12  FOCES  DUE  TO  bOUNDARY  DISP.  FCRA<f  =  l 

D^G  N.y  Ne 

.0(  0.  10  IO0C0  C. 00000000 

3.0C  -.00000009  -.00000011 

6.00  -.OOOooOll  .00000006 

9.00  .COOuGOul  .O0o0ul40 

12. oO  .00000063  .00000432 

15.00  .000OC175  -.oOoio501 

18.00  -.000OU129  -.ooo09095 

-.O00o2575  -.-.,023934 

24. OC  -. 000-06411  .Ool2ol84 

27.00  .00027016  .01130635 

30.00  .00246386  .00113801 
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TABLE  2.IJ  MOMENT  DUE  TC  BOUNDARY  DISP.  FOR  *H  =  1 

DEG  M*  Me 

.00  0.01  OCX  C  0.00000000 

3.00  .00000000  -.OOOOQOOC 

6.00  .00000002  .OOOOOOC 

9.00  .00000004  .00000002 

12.00  -.00000009  .OOOOOOOC 

15.00  -.00000087  -.00000020 

18.00  -.00uu01b6  -.00000060 

21.00  .00001073  .00000167 

24.00  .00007820  .00001820 

27.00  -.00000511  .00001655 

3'  .1  I  -.00277834  -.00055567 

TABLE  2.14  MOMENT  DUE  TO  BOUNDARY  DISP.  FOR  A^= 1 

DEG  Mf  MB 

.1  i  O.OOOOOOI  I  O.OO 

3.00  .OOC       I   .  -.00000000 

6.00  .0OOuC006  .00000002 

9.00  .00000009  .00000006 

12.00  -.00000051  -. 00 o 00006 

15.00  -.  00  0003 16  -tbi.vii.u81 

18.00  -.00000291  -.UOU0U157 

21.0C  .00005065  .00000924 

24.00  .L0C2559O  .00006337 

27.00  -.00045294  -.00003281 

30.1  '  -.01076437  -.L.L224308 
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Diagram  of  the  forces 


-8094  lbs/ft 


-6495.19  lbs/ft 


-6919.9  lbs/ 


■6495.19  lbs/ft 


Dia.  2-1  Dead  load  membrane  force 
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-3237.6   lbs/ft? 


-2767.95   lbs/ft 


•2593.08   lbs/ft 


•2598.03  lbs/ft 


Dia.    2-2     Live   load  membrane   force 


Dia.  2-3  Bending  moment  due  to  AH  =  1  at  boundary 
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Dia.  2-4  Bonding  moment  due  to  0  =  1  at  boundary 


9998.26  lb/ft 


Dia.  3-1  Membrane  forces  due  D.L. 


-15073.54  lb/ft 
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5668.14  lb/ft 


-265S8.85  lb/ft 


753546.23  lb/ft 


Dia.  3-2  Membrane  forces  due  to  Dome  L.  P 


MEMBRANE  AND  RIGOROUS  ANALYSIS   FOR  CONICAL  SHELL  WALL 
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Membrane  forces  due  to  uniform  distributed 
load  and  dome  load 


Fig.  3-1   Section  through  the  conical  shell  wall 


From  the  figure  shown  above  r  =  ySina,   the  shell  is  loaded  with 
uniformed  load  q.   Then  the  normal  forces   P  =  -qSina,  and  tangential 
forces  I- _  =  -q  Cosa  per  unit  area,  distributed  uniformly  with  respect  to 
the  axis  of  the  cone.   The  static  equilibrium  will  require  that  along  any 
hoop  circle  at  a  distance  y  from  the  ape::,  it  will  be 


—  (M1  Cosa)y  Sinad0  =  -y(F  Sinad0Cosa  +  P  SinadB) 
dy   y  y  z 


Which  gives 


/   (P  Tana  +  P  )ydy 


(3.1) 


And 
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,,i  q     Sin  a  ,    .L  q  y2  C 

Nv     =    ~(  rr,,.r     +  Cos  a)  /     ydy      =     — ; +-  (3.2) 

y         y    Cos  a  •  2ycosa      y 

Assume   the  base  of   the   cone    is   fixed,    N'    =  0     at     y  =  L.      Hence 

y 
,    .    .    "' 


2  Cos  a 


2  2 

,  q   L     -  y 

The  hoop  force  N  can  be  derived  from  the  equilibrium  condition  in  a  direc- 
tion perpendicular  to  the  surface,  which  gives  directly 

M'   =  .  p  yTana  =  q  y  Sina  Tana  (3.4) 

For  concentrated  dome  load  P  lbs  per  liner  foot  acting  along  the  top  edge, 
the  static  vertical  equilibrium  equation  is 


M  ^ny  Sin  a  Cos  a  +  R     =     0 

,  -R 

y     "    2rrySinaCosa  C3.3) 


VJhere   R  is   total   loading  due   to   P 

Using  the   equilibrium  condition   in  a  direction  perpendicular  to  the   surface, - 

yields 

,  g  Tan  a  RTana 

Ne    =   5n7?^'ySina  ■=   -TT"  (3-6) 

Using  formulas  (3.3),  (3.4),  (3.5)  and  (3.6)  the  membrane  forces  due  to  dead 
and  dome  load  acting  on  the  shell  can  be  calculated.   For  the  numerical 
values  see  table  3.1  and  3.2. 
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TABLE  3.1  MEMBRANE  FORCES  DUE  TO  DEAD  LCAD 


Y 

65.45 
70.45 
75.45 

8  0  .  4  5 

85.45 

90.45 

95.45 

1 0  u  .  4  5 

105.45 

110.45 

115.45 


-13819.70800000 

-11874.37800000 

-10120.6091      OC 

-8522.68420001 

-7053.24710000 

-5690.9893C 000 

-4419.06770000 

-3223.991400C 

-2094.83130000 

-1022.63450000 

O.OOuOOOOu 


5668.13570000 
6101.14840000 
6534.16100000 
6967.17370000 
7400.18640000 
7833.19900000 
8266.21170000 
8699.22430000 
9132,23 700000 
9565.24970000 
9998.26230000 


TABLE  3.2  MEMBRANE  FORCES  DUE  TO  DOME  LCAD  P 


Y 

f.,  5  .  4  5 

7  0.45 
75.45 
80.45 

8  5.45 
90.45 
95.45 

100.45 
105.45 
111  .45 
115.45 


-26588.85200000 
-24701.77700000 

-23064.815:-  (  01 
-21631.32  7000C 
-2- 365. 5970;  000 
-19239. 80400000 
-18231.95700000 
-17324. 4430000U 
-16502.99000000 
-15755. 91u00<  00 
-15073.54100000 


M9 
7  535  46.23000000 
753546.23000000 
753546.23000000 
753546.23000000 
753546.23000000 
753546.23000000 
753546.  <;3000000 
?53546. 23000000 
753546.23000000 
753546.23000000 
753546.23000000 
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Displacement  from  membrane  theory 

From  Eq.  (1.16),  for  r  =  <»  ,  r  d0  =  dy  ,  r2  =  yTana  ,  »'  ■  Ni  , 
therefore  obtains 

**  -  pw»;  ■*»? +  *.«*  -  ^  -  %*  ^  -*»;> 

I.e. 

For  the  uniform  load  acting  on  the  shell,  the  £q.  (3.3)  and  (3.4)  can  be 
substituted  into  Eq.  (3.7),  then 

A0  «    «T*Ba  [y(y  -V)  -  L2(l  +V)]  (3.8) 

2Ety  Cos  a 

The  horiaontal  displacement  All,  can  be  derived  directly  from  Eq.  (1.7) 

Substitute  Eq.  (3.3)  and  (3.4)  into  Eq.  (3.9),  thus 

AH     .     2|^[ySiBnttoB  +Y£L^Ll]  (3.10) 

Et         -  2yCosa 

Using  the   same   procedure,  &0  and  AH  due   to  dome   load  which  acts  at   top  edge 
of   the  shell   can  be  derived   as: 

,    ,  E  Tan  a  -  1 

**     =     ^HT-rySinaCosa(V-D-     Tana]  (3.11) 

RySlna,  ,; 

AH     =     ■-,„,-■   (Tana  +  V  r )  (3.12) 

-mtt  ySmaCosa 
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Table  3.3  and  3. 4  show  the  numerical  solution  of  the  Eq.  (3.8)  to  (3.12). 


Bend  jr.?,  of  the  bottom  conical  shell  wall 


Co  so. 


Fig.  3-3  Section  through  the  conical  shell  wall  (Showing  the 
relation  of  the  edge  force) 

The  bottom  wall  of  conical  shell  will  be  regarded  as  consisting  of  a 
large  number  of  longitudinal  beams  supported  on  transverse  elastic  rings 
whose  diameter  increases  in  proportion  to  the  distance  from  the  apex  of  the 
cone.   The  bending  analysis  can  be  treated  as  the  problem  of  bending  of  a 
beam  on  a  elastic  foundation.   The  shell  as  shown  in  Fig.  3-3  will  be 
assumed  to  be  under  th->  action  of  an  edge  loading  consisting  of  horizontal 
force  P  and  moment  H  uniformly  distributed  along  the  edge  circle  of  the 
shell.   The  modulus  of  the  foundation  furnished  by  the  hoop  rings,  per  unit 
length  of  circumference  of  the  rings   K  =  -y  ,  where  r  =  ySina  ,   and 

ro 
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TABLE  3.3  DISP.<j  ROTATION  DUE  TC  MEMBRANE  THEORY  (D.L.) 


Y 

65.45 
70.45 
75.45 

8  '.45 
85.45 

9  r.  .  4  5 
q?,45 

10i  .45 
105.45 
110.45 
115.45 


E*AH 
513331  .59000001 
5400X1.47001  00 
566440. 92000000 
590818.96    C     I  I 
631145. 55000C! 0 
665420.690001  ' 
7'   1644.37.  Di  00< 
739816.60.  Oi  OOO 

779937.40 

822006. 78^00000 
866024. 7000000>J 


L*A+ 

-15512.81400000 

-13577.76900000 
-11841.80000000 
-10267.79100000 
-8827.31020000 
-7498.21220000 
-6262.99430000 
-5107.63700000 
-4020.  78170000 
-2993.12430000 
-2^.16.97350000 


TAbLE    3.4    DISP.   ij    ROTATION    QUE    TO    MEMBRANE    THlORY     1 DCME    L.) 


Y 

65.45 

70'. 45 

75.45 

80.45 

85.45 

9i,.45 

95.45 

LOO. 4  5 

105.45 

llu.45 

115.45 


E*AH 
37441834. OOuOOf Oo 
40267633.00001  000 
43093430.   iOOl 
45919229.!  OOC 
48745028..  OO'OI  000 
515  7U827.00OO00O0 
54396626.00000000 
57222425.00001  I 
60048224.00000000 
62674023. oOuouOOj 
66699822 .00 l OC 000 


I  »A4 

-93775. 96900C00 
-92468. 56300CC0 
-91334.43900000 
-90341.29500000 
-89464.36900000 
-88664. 39200C0C 
-87986.14100000 
-67357.39600000 
-o6786. 27200000 
-86270.68400000 
-05797.92200000 
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d  =  t/Cos a  is  the  thickness  of  the  rings  in  the  direction  normal  to  the 

axis  of  the  cone.   Hence  the  modulus  per  unit  length  of  the  longitudinal 

beams  will  be  K  -  blC  where   b  -   b,y   is  the  width  of  the  beams  increasing 
0-0 

linearly  with  the  distance  from  the  appex.   Thus 


K  =  bnyK 


bQEt 


y  Sin  a  Cos  a 


The  flexural  rigidity  of  the  beam  will  be 

byt3E 

D  .  EI  =  D  b  y  = H-2 1- 

u      12(1  -  V")CosJa 

Putting  these  values  of  1C  and  EI  into  the  differential  equation  of  bending 

2  2 

— 5-(EI— |)     +  KZ      =     0  (3.13) 

dy  dy^ 


,2  b  vt3E  .2,  b  Et 

i °1 LZ  +  2 =     o  (3.14) 

dy     12(1   -V   )Cos  a  dy  ySinaCosa 


After  rearranging  and  differencing  of  equation    (3.14),   yields 


y<L|  +  a^  +  ?Cot2a     .     0  (3.15) 

dy  dy  t   y 


The   finite  difference  equations  will   therefore   be 


Z.,,   -  4Z.  .,    +  6Z.  -  4Z.    ,    +  Z.    „            Z.    „  -   2Z.    ,    +  2Z.    ,    -   Z.    . 

f_j+2            i+l            *  1-1           1-2.         .    i+2             i-H             i-1          i-2 

I  )+2(-                                 -                              ) 

h  2  h 


12  Z   (1   -V2)  , 

+  2* Cot  a     =     0  (3.16) 

t  y. 
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Collecting  similar  terms  of  Z's 


z         Zi       L+!        A        '"yU     ,    -    f6yi     ,    12Cot2q(l-A 


_4y(        2  y-      i 

+  Z.  „  (— 7-1 ^)   +  Z.  „<4  +  4)      =     0  (3.17) 


Let 


i+ls  .4  .3'  H2  \4       .3' 

h  h  h  h 


4  j  1 


h  h 


L     4yi    -    r 
7     7     "    '« 


6yi    .   12Cot2q(l   -  V?) 

v.*  .2  =        1 

h  y.t 


•*y.       2 


4         3  i 

h         h 


Then  Eq.    (3.17)   becomes 


A.Z.    ,   +  B.Z.    ,    +  C.Z..+  D.Z,    ,    +  E.Z,    -      =     0 
i    1-2  1    1-1  i    i  i    i+1  i    H-2 


In  this  problem  h  =  5'  ,   total  there  are  eleven  points.   Therefore  there 
are  eleven  equations  with  fifteen  unknowns.   These  eleven  equations  are: 


A,Z,    +  B,Z-   +  C.Z.,   +  D,Z,    +  E-Z.      *     0 
31  32  ]   J  34  35 
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A.Z.    +B.Z,    +  C,  Z,    +D.Z=   +E   Z,      =     0 
42  43  4  4  45  6 


(3.18) 


A13Z11   +  B13Z12   +  C13Z13   4  D13Z14  +  E13Z15     =     ° 


Assuming  zero  displacement  and  rotation  at  bottom  edge,    (point   3) 


z3     =    0 


z,  +  z4    =    0 


(3,19) 


Assuming  unit  horizontal  displacement  and  zero  rotation  at  top  edge  (point  13), 
gets 


Z    =2 

13 


7    -!-  7 
12  '  Z14 


(3.20) 


Fig.  3-4  The  division  of  the  conical  shell  wall  in  force  analysis 
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By  moans  of  above  assumption  and  combining  Eq.  (3. IS).  (3.19)  and 
(3.20),  there  are  a  total  of  fifteen  equations  with  fifteen  unlcnowns.  Hence 
these  equations  can  be  solved  and  get  the  displacement  of  each  point  corres- 
ponding to  the  unit  horizontal  displacement  at  top  edge  (point  13).  Further- 
more, to  use  these  value,  calculate  the  stiffness  of  the  beam  due  to  this 
displacement. 

Uith  the  same  concept,  assuming  bottom  edge  fixed  and  zero  displacement 
and  unit  rotation  at  top  edge,  it  Is  seen  that 


z3  =  0 


z2  +  z4  =  0 


Z13  "  ° 


hi   +  Z14  =  10 


(3.21) 


Eq.  (3.21)  together  with  Eq.  (3. IS)  will  provide  another  set  of  solutions. 
Then  the  stiffness  of  the  beam  due  to  rotation  can  be  solved.  Using  the 
same  procedure  the  stiffness  at  bottom  edge  can  be  determind. 

In  solving  the  above  equations,  use  the  same  method  as  before.   Write 
the  equation  in  matrix  form 
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A3   B3   C3   D3  E3   0 


0     A,    B.    C,    D,    E, 
4     4     4     4     4 


0     0     A.   B,   C.  D. 
5      5      5      5 


0  0  0  0  0  0 
0  0  10  0  0 
0-10100 
0  0  0  0  0  0 
0     0     0     0     0     0 


G2        G3 


0  0  0  0 
0  0  0  0 
0        0        0        0 


B13   C13  D13  E13 
0       0       0       0 

0       0       0       0 

0       10       0 

-10        10 


zll 

z12 

z 

Z14 

UsJ 

0 

~  0" 

'0' 

0 

0 

D 

0 

0 

0 

0 

0 

• 

or' 

for' 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

10 

2 

0 

0 

0 

or 

.0, 

.0, 

Let  the  coefficient  matrix  be  A,  the  Matrix  of  constant  term  of  the  equation 
be  G,  then 


(z) .  mX 


The  A,  A-  and  the  solution  of  those  equations  see  table  (3.5),  (3.6),  (3.7), 
(3.8)  and  (3.9). 

Forces  due  to  cde;e  effect 

From  above  analysis  there  are  already  obtained  the  displacement  due  to 
edge  effect.   Kou  use  this  value  of  displacement  to  calculate  the  bending 
moments  and  forces  both  at  meridional  and  hoop  directions.   From  general 
theory,  the  moment  is  equal  the  second  direvative  of  displacement  in  Z  direc- 
tion with  respect  to  y,  therefore  the  moment  at  meridinal  direction  can  be 
obtained 


TAbLE    3.5    COEFFICIENT    MATRIX    CF    AX=G 
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.09672o 

.11272' 

.000 0U€ 

.000000 

.10472. 

.120720 

.000000 

.000000 

.112721 

.126720 

.000.     ! 

.000000 

•12^72o 

.136720 

i00«    iOi 

.ooc 

.12S72- 
.144720 

.000000 
.i  iQi  K 
.136721 
.152720 
.000000 
. C 0 000 c 


-.402880 

.1  ooooc 

.    I 

0.000000 

-.43488;.' 

.       OOC 
Q . 00  00 ( 
-.46688C 
0.00000' 
0.00 

i..  uOOOOO 
-.498680 
C.OC 
0.  OOC  ('0  I 

.  '    ;       0 1 

-.53086'.: 
0.00 
.     C     I 

-.562880 
O.uOOOOo 
D. OOOOC 
0.000000 


.760329 
o  •  ,  i  ■  0  0  0  0 
0 . 0 00  u(j^ 
0.0  I U  0  0  0 

.798960 
(  .  ■      OOC i 
0« 
'  .         1 0  0  0 

.838833 
o. 0  0 0000 
o  .  i  OC  OOC 
u .000000 

.879716 
0.000000 
O.OOC  0<  0 
0.000000 

.921432 

u. 

OOo 
0.000000 

.963642 
0. jOOOOO 
0.  [ 

o.uOoooo 


-.434880 
0.000000 

0.000000 
0.000000 
-.466880 
0.000000 
0.000000 
0.000000 
-.498880 
0.000000 
0.000000 
0.000000 
-.530880 
0.000000 
0.000000 
0.000000 
-.562880 
.00 3000 
0.000000 
0.000000 
-.594880 
0.000000 
0.000000 
0.000000 


.14472 

.160720 

.000000 

.000000 

.152720 

.16o72l 

•  01  OOG 
.160720 

.176720 

..'-"<;■.■'.( 

.000000 

.168720 

.184720 

.001..  )0 

.0001  00 

.176720 

.192720 

.00( 

.000000 

-1.000000 

.000000 

.000000 

.000000 

1.0  0 ( . 

-l.OOC 


-.594880 

.001 

.1  OOOOO 
O.OOOOOCi 
-.626880 
0.000000 

o.t.ooooo 
.  ooooo 

-.658880 
.  OOOOO 

.  ooooo 

O.OOl 
-.690880 
I  .00000 
0.00000 j 

.      .      iO 
-.  ^22660 
0.000000 
.  .000000 
0.01    lOi 
0.00000 

0.  ooooo 
o.oooooo 

0.0 ooooo 
! .000000 
O.OOOOOO 


1.006839 

O.00  Ol 
0 • uOO  0  0  0 
1.050333 
0.  000000 
u.GOQOQO 
o.OOOOOO 
1.094255 
0 .  . .  OOOI 
0.  )0000 
o.<  OOOOO 

1.138545 

0.01      00  C 

0  .  (  I  '...  0  0  o 

O.uOUOOl 

1.183158 

u.uOOuOo 

0.00 OGOu 

0.000000 

I..        Q< 

0.00 

O.OOUC 

0.      (OOOI 

0. 

1. jOOOOG 
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- 

.626880 

0 

.000000 

0 

.000000 

0 

.000000 

- 

.658880 

0 

.000000 

0 

.000000 

0 

.000000 

- 

.690880 

0 

.000000 

0 

.000000 

0 

.000000 

- 

.722880 

0, 

.000000 

0, 

.000000 

0, 

.000000 

-< 

.7548B0 

0. 

.000000 

0, 

,000000 

0, 

.000000 

0.000000 

0, 

,000000 

0.000000 

0. 

,000000 

0. 

,000000 

0. 

,000000 

.000  0.00                                  0.0001                              0.000000 

.000000  0.0001                              O.OOl                                 0.000000 

.000000  0.000000                      i  .       i  01                          0.000000 

1.000000  .   .001  :  I                            0.1      I  0                             0.000000 

•  OOOOOi  0.00000                      o.soooop,                   0.000000 

.000001  .      '0000                       0.00                                  0.000000 
.000    - 

TABLE  3.6  INVERSION  OF  MATRIX  A 

.10339123E+02  -0.27975251 E+01  -0. 3016985CE+C 1 

.131  62459E+G2  0 . 742O6158E+01           0. 26669520E+0 1 

.3030497b  -C. 44017222  -0.45849873 

.27854976  -0. 11626985  -0. 25 19933 lfc-Ol 

.00000000  0.46548200E-02           0. 23346486E-02 

•  uOGui  -C. 82443642                 0.58797084 

.16765051E+01  0.12518429E+01           0.57518303 

.15216664  -0.22554u37E-0l  -0. 60050102E-0 1 

•45649217E-CI  -0 . 22363885E-01  -0. 58558356E-02 

.OOOOi  D.10816901E'-02  -0.28092077E-03 

B-.000000  0.10000000E  +  01 

.0      lOOOl  ..  OC     0                                                0.00000000 

.0000001  0.000000                                             0.00000000 

.000001  O^OOOOIX                                              0.00000000 

1   .00000000  0.00000000 

.OOOOOUOU  0.175D6360                                        0.58797084 

•16765u51E+01  0.12518429E+01                             u. 57518303 
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.1 5216664 

.45649217E-M 

.001  01 

.0  10    '  i  i. 

•12518428E+01 

.87090266 

.10026390 

. 

.57518301 

.22291010E+01 

.46048117E-01 

.2.  I000OOE-  8 
. 

.1521666? 
.3515816C E+01 
.3: 934961 

•  ■ 

•  000  00  0  00 
•22554U35E-01 
.23523037E+01 
.11719418E+J1 
.( 

. 

•60050099E-01 
.  10763079E+01 


-! .22554037E-01 
-: .22363885E-01 

C  .108169C1E-02 

■   .13109301 

i  .29366761E  +  01 

C. 210^6576 
-C.6643u234E-cl 

0.4Q452676E-02 

0.6G233167E-01 

G.19864946E+GI 

C.9bl25517 
-0.895O5870E-01 

C.75476550E-02 

0.15934886E-01 

0.87090267 

C .23523038E+01 

0.48937200E-02 

0.68139390E-02 
-0.23618594€«»ef2 

C.210J6578 

0.36065352E+01 

0  .  37227  u 10 
-0.8300815QE-O2 
-L.62884465E-02 
-0.52888184E-01 

C.2463u225E+0l 


-0.60050102E-01 

-0.58558356E-02 
-0.28092C77E-03 

0.21337940 

0.19864947E+01 
-0.52868221E-01 
-0.21899456E-01 
-0.40911268E-02 

0.26217861E-01 

0.33770157E+01 

0.24497700 
-0.40859960E-01 
-0. 13719374E-01 
-0.31993933E-01 

0.22291011E+01 

0.10763078E+01 
-0.36887910E-01 
-0.27530351E-01 
-0.33486908E-C1 

0.98125516 

0.24630225E+01 

0.44937280E-01 
-0.33303301E-01 
-0.20153025E-01 

0.24497700 

0.36726315E+01 
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.Z508B310E+01  I . 1 1475138E+01  0.27923684 

.100OOO00E-09  -- .5158O633E-01  -0.93392267E-03 

.00000000  -6.47803859E-02  -0.8550182 1E-02 

•45649216E-01  -C. 10026388  -0. 46C48 1 1 7E-0 1 

.30934961  L.11719418E+01  0.25088310E+01 

•35575475E+01  0 . 221 85361 E+C 1  0.70470680 

.200000005-07  -0.13017348  0.11735843 

-(  .23419459E-02  -0.22375902E-02 

•22363884E-ul  -0 .6643u230E-L  1  -0. 89505880E-01 

•4893670OE-02           0.37227011  0. 11475138E+01 

•22185381E+J1           C .283bl276t+01  0. 1 1989235E+01 

.00000000  -0.22146511  0.36512378 

jOOO  -0.61322310E-03  -0.78079200E-04 

•  5855836UE-J2  -v.  .  2  1899460E-C1  -0.  40859970E-01 

.36887930E-U1           C  .44937296E-01  0.27923690 

.7047068C               C.11989233E  +  01  0. 1 2918149E+01 

-L. 23862410  0.72195348 

.OOOOi                                                 .00000000  '."0000000 

.001  OOOOo                 .  0(  0.00000000 

00000  •'         IJ(  -0.  10000000E-06 
.    0000001  0.00000000 

O.OOOuuUOC  .  *  .  I 

00001  -0.61322310E-03  -0. 78079200E-04 
•5855836uE-u2  -0.2 1899460E-01  -0. 40859970E-0 1 
,3688793uE-ul           0.4493  7296E-01  0. 27923690 

.7  47                  0.11989233E+01  C. 12918149E+01 
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.  i ■■■:.!  0.72195348 

.0000'                  -  .25546262E-02  0. 14531 162E-02 

.24394844E-01          -0.10700803  -0.23123572 

.28734082               0 . 32 1 17190E-02  0. 1 0889186E+01 

.3369279?E+ol           0.659O7381E+O1  0. 88061370E+01 

.51886751t+oI           0.229U3U86E+O1  -0.93819192 


TABLE  3.7  CONSTANT  MATRIX  G 


5| 

o.O 

10.0 

O.o 

O.O 

.0 

o.O 

O.o 

O.O 

.0 

0.0 

O.o 

G» 

. 

o.oc 

2.o 

.0 

CO 

0.0 

.0 

u. 

O.o 

G, 

. 

0.0 

. 

.o 

.0 

o.O 

Co 

o.O 

. 

10>C 

O.o 

G1 

0.0 

o.O 

CO 

U  •  u 

.0 

o.O 

0.0 

~  »0 

.0 

o.O 

2.o 

0.0 
0.0 


o.O        0.0 

o.o        0.0 


0.0 

0.0        0.0 


o.o        0.0 
o.O        0.0 


TABLE  3.8  W  VECTOR  FOR  TCP  EDGE  FIXED  A^>  = 1  CRA"=1  AT  BOTTOM  EDGE 

Y  A<f  =1                                              *H"1 

55.45  -27.975251  C          -6.03397000 

60.45  -8.24436420           1.17594170 
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65.45  .  2,00000000 

70.45  1.75563'  1.17594170 

75.45  1.31093010  .42675880 

8o.45  .6*-233167  .05243572 

85.45  .15934886  -.06398787 

90.45  -.02361859  -.06697382 

95.45  -.00288447  -.u40Ju605 

lou.45  -.04760386  -.01710036 

103.45  -.02341946  -.00447516 

110.45  -.00613223  -.00015616 

115.45  .    00000  .      i  00 

120.45  -.00613223  -.000L5616 

125.45  -..2554628  .00290623 

T4BLF  3.9  W  VECTOR  FOR  BOTTOM  EDGE  FIXED   A+=l  0RAH=1  AT  TOP  EDGE 

A+  =  l  4H»1 

55.45  .04654820  .00466930 

6o.45  .01081690  -.00056184 

65.45  U.0bO;OOOOO  U.uuuooOOO 

70.45  .01081690  -.00056184 

75.45  .04045268  -.00818225 

80.45  .07547655  -.o2743875 

83.45  .06813939  -.05506070 

«o.45  -.08300815  -.06660660 

95.45  -.51580633  -.00186785 
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I  C.45 

I I  5.45 
110.45 
115.45 
120.45 
125.45 


-1.3,:  173480 
-2.214651 10 
-2.38624100 

0  .  0  C  0  •      I 
7.61375900 
22. 9^308600 


.23471686 

.73024756 

1.44390700 

I  .   ooooo 

1.44390700 
-1 .07636380 
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H27 
M   =  -  Dn  ~  (3.21) 

y        °dy2 


Where  D   is  flexural  rigidity  of  the  sh3ll 


Et3 

%  =  — H — —  °-n) 

0     12(1  -  V  )CosJa 


The  hoop  bending  moment  will  be 


ti   =  VM  0.22) 

o     y 


The  shearing  force  Q  can  also  be  expressed  in  terms  of  the  meridional 
bending  moment  M  by  a  consideration  of  equilibrium  of  the  meridional  beam 

y 

(see  Fig.  3-3). 

3     2 
dyMv         d  Z   d  Z.  ,.,  ,,, 

Q   =  -  -T~   =  Dn(y— 1+— 9)  (3.24) 

y        dy       °   c!y3   dy2 

The  meridional  force  K  can  be  obtained  as  a  component  of  the  shearing 

y 

force  Q  : 

y 

3     2, 

N   =  Q  Tana  =  D  Tana  ( — T  +  — f)  (3.25) 

y     y  °      dyJ  dy* 


The  hoop  force  N  will  be  proportional  to  the  deflection  Z  and  according  to 
Hook's  law,  its  values  per  unit  length  of  the  generator  will  be 

Ft 

Kn  = ^— Z,  (3.26) 

0     y.  Tana  1 

Again  write  the  equation  (3.21),  (3.23),  (.24)  and  (3.25)  in  finite  dif- 
ference 
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Z.   -  2Z,  +  Z. 

y  °  b2 


Z    -  2Z  +  Z 


Z  _  -  2Z    +  2Z.    -  2.  ,    Z    +  2Z.  +  Z.  , 

Q   =  V" W  ■  M ^  +  -^ ^— ^)    (3.27) 

y      °  2  li  h 


z,  n   -   2zi.u  +  2Z-  i  -  z-  •>   z-  i  "  2Z*  +  Zi.i.i 
,  i-l-2     i+I     i-I    i-2    l-l     l    1+1  . 

N   =  D.TanaC 5 + S ) 

y       °  2h3  h2 


Using  the  data  of  Z  value  obtained  from  former  calculation,  the  M  ,  H  ,  Q  , 
N  ,  N  can  be  obtained  as  shown  in  table  (3.10)  to  (3.13). 
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TAfciL       3. It    FORCES    DUE    TO    EDGE    EFFECT     FOR    TCP    EDGl    FIXED 
AND     AH=u,Af=l    AT    BOTTOM    LUCE 


0 


6 


■ 


d 


I 


65.45  1.95377680    3.3841   397C  CO    OOOOOC       •■    I010C94       .00002019 

70.45  .835^0296  1.4467867C  .00959184   .02B99987   .00579997 

7*. 46  .16174812   .28015588  .0'  66F756   .02167972   .00433594 

80.45  -.07441295  -.12688698  .0-288176   .01001572   .00200314 

85.45  -...9731782  -.16855936  .0u071777   .00272393   .00054479 

9  0.45  -.05916194  -.10247146  -.0-010051  -.00028784  -.00005757 

95.45  -.02314657  -.04009102  -.00025358  -.00093418  -.00018684 

100.45  -.00344401  -.00596519  -.Oo018317  -.00068594  -.00013719 

1-5.45  .00343091   .0u594251  -. 0^.008548  -.0OU28456  -.00005691 

110. '5  .00373356   .0^646671  -.0*002137  0.      000  0.00000000 

115.45  .00192227   .00332947  O.OoOOOOOu   .00010094   .00002019 


TABLE  3.11  FORCES  DUE  TC  EDGE  EFFECT  FOR  TCP  EDGE  FIXED 

AND  AU=l»Af=o  AT  BOTTOM  LDGE 

Y        J            Hq                          H9  M+                      Mj 

65.45   .67395451  1.17078720:   .01176167  .03292389   .00658478 

70.45   .05940799   .10289763   .0-642471  .01935935   .00387187 

75.45  -.10623272  -.18400041   .00217707  .00702730   .00140546 

80.45  -.09261803  -.16041909   .00025087  .0uu86570   .00017314 

85.45  -.04565829  -.07908246  -.00028823  -.00105071  -.00021014 

91  .45  -.012913'  5  -.02236606  -.Oo02850>J  -.0ol09986  -.00021997 

95.45   .00150582   .00260816  -.00016253  -.00u66o89  -.00013218 

i   .45   .00490867   .00850207  -.0.006552  -.00027891  -.■Q0005578 
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105.  45  .0039499  .(  684143  -.  0  01633  -.  007109  -.00001422 
110.45  .00224562  .00388952  -.00000054  u.  0GOOOL  .  ...0000 
115.45   .00112403   .00194688  O.Ov*  )00     .      257   .  I  00051 

TABU  3.1?  FCRCti,  DOE  TC  EDGE  EFFECT  FCR  OCTTCM  EDGE  FIXED 

AND   4H=v.^  =  l  AT  TCP  EDGE 

Y  N^  Nj  *U  Mj 
65.'  5  -.00265773  -.00460333  O.i  -.  i-302420  -.00860484 
70.45  -.00132525  -.00229540  .0u00591u  -.U4284615  -.00856923 
75.45  -.00342186  -.00592683  .00020637  -.4235833  -.00847167 
80.45  -.01032255  -.01787917  .00036111  -.04178181  -.00835636 
85.45  -.02029657  -.03515467  .'  10693  -.04190258  -.00838052 
90.45  -.  :  . 2  246  -.04191968  -..  -.04439057  -.00887811 
95.45  -.00  035446  -.00061398  -.0u207998  -.5151472  -.01030294 

1  i  .45  .08558652  .14824015  -.00498794  -.06445102  -.01289032 

105.45  .26565260  .46ul2367  -.0,;B08304  -.07947863  -.01589577 

110.45  .527?5795  .91375691  -.00831567  -.00230331  -.01646066 

115.45  . 739/4176  1.28126990  0.0^1.0000*.  -.04302420  -.00860484 

TABLE  3.13  FCRCti  DUE  10  EDGE  EFFECT  FCR  BOTTOM  EDGE  FIXED 
AMD   AH=li*t  =  C  AT  TCP  EDGE 

Y  N*  Ng  M>         M« 
65.45  -.00095305  -.00165074   .     .   !;   -.2376767  -.00475353 
70.45  -.00063968  -.0110796  -.(    I   !>  1    -.  2377691  -.00475538 
75.45  -.00024427  -.0   423  8  -.'   04174  -.02390235  -.00478047 
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157743  .0i'273218  -.00013126  -.02421933  -.00484387 

85.45   .00620787  .01075235  -.   024802  -.02467400  -.00493480 

.45   .01371832  .02376082  -.00028344  -.02486405  -.00497281 

95.45   .02007407  .03476929  -.00000753  -.02379841  -.00475968 

100.45   .01339599  .02320253  .00089938  -.01990407  -.00398081 

105.45  -.02818109  -."4881106  .0  266546  -.01174731  -.00234946 

110.45  -.13168:165  -.22807751  .0  .'503178  i..'.t"".vi'  0.00000000 

115.45  -.30976045  -.53652068  .Oi 666783   .00915366  .00183073 


75 


Dia.  3-3  Top  edge  fixed,  and  AH  =  0,  t><t   -   1 
at  bottom  edge 


The  Ola.  on  above  and  the  following  pages  are  just  for  the  purpose  of 
showing  the  pattern  of  forces.   Therefore  no  scale  and  dimension  are 
presented. 
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Dia.  3-5  Top  edge  fixed  and  Ali  =  1,  40  -  0  at 
bottom  edge 
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-.Dia.  3-6  Botton  edge  fixed,  and  AH  =  0,  A0  =  1 
at  top  edge 


7  a 


Dia.    3-7      Bottom  edge   fixed,    and  AH  =  1,  A0  =  0 
at  top  edge 
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THE  CALCULATION  OF  THE  FLEXIBILITY  OF  THE 
DOME  AND  THE  CONICAL  SHELL  WALL 


In  the  previous  section  the  stiffness  of  the  dome  and  the  conical  shell 
wall  had  already  obtained.  The  next  step  is  to  apply  the  reciprocal  lau  to 
calculate  the  flexibility. 

At  first  the  dome  is  considered.  Based  on  statically  equilibrium  con- 
dition, for  unit  horizontal  displacement  and  zero  rotation  at  edge  of  dome, 
we  have 


FUH11  +F12H21   =  1 


F21HH  +F22H21   "  ° 


(4.1) 


Khere  H  represent  stiffness,  F  is  flexibility  for  correspond  ins  relation, 
see  Fig.  4-1. 


Fig.  4-1   The  relation  of  stiffness  &  flexibility 


Then  assume  there  is  o  unit  rotation  and  a  zero  displacement  at  the 
edge,  hence 
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Fll  "l2     +    F12  H22      "     ° 


F2l"l2     ^    F221!22      =     l 


(4.2) 


Combining  Eq.  (4.1)  and  (4.2),  yields 


Hence 


Fn     Fi2\   /Hn     Hi2\       />       ° 


VF21        F22/     \H21        W  \° 


Fn     Fi2\       [hi     "12V1 


'F21   F2?/    \n21   H22, 


(A. 3) 


(4.4) 


Using  the  same  idea  and  procedure  the  flexibility  matrix  of  the  bottom 
conical  wall  can  be  obtained.   For  AH  =  1,A$  =  0  at  top  and  let  the  bot- 
tom of  the  wall  fixed,  then 


F33H33  +  F34H43  +  F35H53  +  F36H63   "  J 


43  33    44  43    45  53    46  63 


F53K33  +  F54H43  +  F55H53  +  F56H63   =  ° 


F,,H,,  4  Fr.H,,  +  F,-H-„  +  F,,H,,   ■=  0 
63  33    64  43    65  53    66  63 


(4.5) 


For  AH  =  0,  A0  =  1  at  top  and  still  let  bottom  edge  fixed,  obtains 


F,,H,,  +  F,.H.  .  +  F,.H.,  +  F,,H,,   =  0 
33  34    34  44    35  d4    36  64 


F  H   +F..H..  +  F.-H..  +  F,,H,, 
43  34    44  44    45  54    36  64 


(4.6) 


F,-,H   +  F  II   +  F-.H.,  +  F.,H,, 
53  34    54  44    55  54    56  64 


F  H   4  F  1-1   +  F  H   +  P  H    -  0 
63  34    64  44    65  54   "66  64 


For  AH  =  0,  A0  -  0     at  top  and  AH  -  1,  .40  =  0     at  bottom  of  the  conical 
shell  wall,   hence 


F33H35  +  V«  +  F35H55  +  F36H65     "     ° 


FH       +FH       +  F     H       +FH  =0 

43   35  44  45         45   55  36   65 


F„    H_   +  F     H,  „   +  F„„H««   +  Fr    K  c      =     1 
53   35  54  45  55   55  56  65 


FH       +FH       +FH       +FH  =0 

63   35  64  45  65   55  66   65 


(4.7) 


For  AH  =  0,  A.0  =  0  at  top  edge  and  AH  =  0,  A0  =  1  at  bottom  edge  of  the 
conical  shell  wall,  therefore 


F,,H,,  +  F,,H.,  +  F..H.,  +  F,,H    =  0 

33  36    34  46  35  56    36  66 

F  H  +  F  H  +  F  H  +  F  H    =0 

43  36    44  46  "45  56    46  66 

(4.8) 

F   H  +  F   H  +  F   H  +  F   H    =0 

"  53  36    54  46  55  56    56  66 

F  H  +  F  II  +  F  H  +  F  H    -  1 

63  36    64  46  65  56    66  66 


From  Eq.  (4.5),  (4.6),  (4.7)  and  (4,8)  there  are  sixteen  simultaneous 
equations  with  sixteen  unknowns.   Therefore  it  can  be  solved  easily.   Write 
it  in  matrix  form: 


82 


F  F  F     ' 

*0i      '35      r36 


33        3 


P.        F         F         F 

43       44       45       46 


F  F  F  F 

53        54        55        56 


F  F  F  F 

63        64        65        6! 


rH33      H34     H35      H36l 


!.,      H  H  \\ 

43       44       45        46 


H„      H,.      H„      E,r 
53        54       5:>        36 


-!  H  H  H 

63        64        65        66j 


>  °S 


I  0        0       0 

0  10        0 

0  0        10 

0  0        0        1 


Hence 


F  F         F  F 

33        34        35       36 


F  F  F  F 

43        44       45       46 


F  F  F  F 

53        54        55        56 


F  F  F         F 

63        64       65        66 


Si  H  I!  HI   "1 

"33      M34       35        36 


"43      H44      H45      "46 


'53      H54      "55      "56 


H63      "64      H65      H66 


(4.9) 


The  value  of   H's  are  already  calculated    in  previous   section.      For  the 
solution  of  Eq.    (4.4)   and    (4.9)   see   table    (4.1). 

TABLE  4.1      VALUES  OF  FLEXIBILITY 


11 


l22 


33 


36 


64 


53 


56 


63 


66 


0.24260895E+05 
-0. 13341623E+03 
-0.92809630 
-0.23932046E+02 

0.51200690 
-0.10017272E+02 

0.92809668 
-0.64433998E+01 
-0.32085288 

0.22039434E+01 


F12  0.48091368E+03 


F,,        -0.41100000E-02 
34 


F, ,        -0.20574000E-01 
44 


F_.  0.30379565E+O2 

54 


F,,         -0.10510547E+02 
64 


F,,        -0.6261863E+04 


FJ5  0.17373000E-02 


F,  =        -0.15800000E-03 
45 


F  -0.17375700E-02 


F65  0.34182137 
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COMBINATION  ANALYSIS  OF  THE  DO'E- 
RIHG-CONICAX  SHELL  WALL 


Due  to  membrane  analysis  both  the  dome  and  the  shell  wall  Induced  a 
great  deal  of  displacement  and  rotation  along  the  edge.   In  order  to  res- 
traint these  displacements  and  rotations,  a  ring  will  be  provided.  This 
ring  acts  as  a  circular  tension  tie.  And  for  the  purpose  of  reducing  tho 
moment  at  edge  of  the  shell,  it  may  logically  be  prestressed. 

Let  the  ring  be  in  the  rectangular  shape.  There  are  three  forces 

ccting  en  it,  M  ,  H,  and  M  as  shown  in  Fig.  5-1.   Corresponding  to  those 
a        x 

forces,  the  horizontal  displacement  H  and  the  rotation  0   are  therefore 
induced.   The  relation  of  it  are  easily  derived.   The  result  are 


H  =  j-H 

Ebd 


H  = 


I2r2y, 
Ebd2 


(5.1) 


12r 

Ebd 


3  a 


Fig.  5-1  The  relation  of  the  forces  which  acting  on  the  ring 


84 


The  relation  among  dome  ring  and  shell  wall  arc  shown  in  Pig.  5-2. 

*R   KK   *D 

620  \o    610      V. 


F  ;* 


Fig.  5-2  The  relation  among  the  dorie-ring-sholl  wall  and  the  sign 
convention 
Fig.  5-2a  and  b  illustrate  the  system  in  which  all  stress  resultants 

for  the  dome  and  shell  wall  are  determined  by  the  membrane  theory  where  N1  . 

ad' 

N  ,  and  K   are  meridian  forces  at  edges  of  dome  and  wall.   H   is  the  rin» 
aw      aw  f 
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force  due  to  prestresslng.   In  this  problem  the  ring  and  dome  supported  on 
the  conical  shell  wall  will  settle  uniformly  every  where  due  to  uniformly 
distributed  load.   Therefore  it  can  be  stated  that  the  vertical  settlement 
has  no  effect  on  the  analysis  of  the  relation  between  dome  and  ring.  The 
effect  of  the  vertical  deflection  of  the  bottom  conical  shell  wall  is 
already  induced  into  the  horizontal  displacement  63  and  6j  (Fig.  5-2c). 

Now  there  will  be  ten  components  due  to  translation  and  rotation  of: 

the  dome,  the  top  and  the  bottom  of  the  ring,  the  top  of  the  shall  wall  and 

D    D    S    R    R    R    w    w    w 
the  bottom  of  shell  wall  (6^,  6^,  6^,  6^,  5^,  6^,  6^,  6^,  6^  and 

6W  respectively).  Fig.  5-2c  illustrates  this  system  and  shows  that,  in 

60 
addition  to  all  the  previous  displacements,  there  may  be  an  additional  ring 

rotation  for  the  center  line  of  the  wall  does  not  intersect  the  centroid  of 

the  ring  cross  section 

aR       6r2e  „t  R   _  !2rV  .T 

2°  =  "  Ebd2  'aD  40     Ebd3   aw 

and  ,  i  2  i 

5    =  f Cos  a  h t-) =  (Cos  a  +  -r-)— — 

10  d2  Ebd  d  tbd 

Where  a  is  the  angle  between  the  radius  of  curvature  and  the  axis  of  revolu- 
tion at  the  edge  of  the  dome. 

12y_c  r2M,T 

,R      ,...     .°  w, aw 

&         ~     (Since  +  ,— ) 

30  c,2    Ebd 

For  6  ,  620,  6j0,  o"0>  6^0  and  6gQ  arc  already  calculated  in  the  section  of 
membrane  analysis.   Then 

610   "   610  +  610 
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620  "  620  'r  °20 


D     R 
630  "  630  +  630 


D     R 
640  =  640  +  640 


In  order  to  restrain  these  displacement  due  to  membrane  theory,  there 
are  four  corrections:  a  forces  X  and  a  moment  X„,  which  correspond  to  the 
required  dome  -  ring  values;  and  a  force  X  ,  and  a  moment  X  ,  which  corres- 
pond to  the  ring—  wall  values;  and  a  force  X  ,  and  a  moment  X  ,  which 
correspond  to  the  wall  and  the  basement.  The  correction  displacements  due 
to  those  forces  as  shown  in  Fig.  5-2d.  Consider  first  the  horizontal  dis- 
placement at  the  junction  of  the  dome  and  ring,  due  to  the  force  X  ,  obtains 

2r2 
*il   EM''  1 

2 

612  "  621  +  62\   =  (F2!  --^2)X1 

Ebd 


From  X  force,  yields 


D  R  6r2 

612      =      612+512      =      (F12-— 2>X2 

Lba 


2 

622     =     622   •'    622     =     (F22  +  I~3>X2 

Eod 
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And   from  the  ring-wall   forces  X     and  X.    oome  the  displacements: 

6        -     6R  A3       12r2y0(d/2)X3 

13  °13  "    Ebd    "  Ebd 

2  2  2 

r  X^       3r  X3  2r  X3 

TbTi ' 


and 


14 


Ebd 

Ebd 

c  • 

"       Ebd3      * 

6r2X, 

4 

Ebd 

The  rotation  of  the  ring  at  this  junction  due  to  X,    is 

12r2y0  12r2d/2  6r2 

23  rhH3    X3         '  p.  ,3    X3       =     *  _,  ,2*3 

Ebd  Eod  Ebd 


and   the  rotation  due   to  X,    will   be 
ft 


a      -      ^x 


In  a  similar  manner  the  displacements  at  the  ring-wall  junction  can  be 
derived. 


631 

= 

2r2X 
1 

Ebd 

632 

= 

12r2d/2 

_   ....              -V 

Ebd3    "2 

12r  . 
Ebd 

S33 

= 

R           w 
633    H    633      = 

2                 2 
rX3        6ry0X3 

Ebd              ,.:bd2 

= 

,   2 

F    x    -  2i-x 

33   3        EbdA3 

+  F,,X„ 


and 


,I!  ,w  6r2 

634     =      634   +  634     "     —: 2  X4    *  F3A 
Eba 


41 


Ebd3 

6r2 
"     Ebd2"1 

-I^X 

642 


R  W  6r2 

643     "     643  +  643     "     ^TX3    +  Wl 


2 
R  W  1  ?  *• 

44  44  44  „  ,3    4  44  4 

Ebd 


645     =     V5  ;  654     "     F54X4 


653      =      Vs  !  646      =     V6 


664  =  F64X4  !  636  =  F36X6 


663  =  F63X3  !  635  =  F35X5 


Now  there  will  be  six  simultaneous  compatibility  equations  to  solve  for 
the  six  corrections  J^,  X,,,  Xj,  X  ,  X  and  X  .   Those  six  equations  are: 


610  +  611  +  6!2  +  613  +  614 


620   +   621   +   622   +   623   +   624 


*,„  +  $„     +  6  „   +  6       4  6       4  6        1-6 
30  31  32  33  34  35  36 


8,  „  +  6.      4  6       +6       +6       +6       +6 
40  41  42  43  44  45  46 


(5.2) 


650  +  V,   +  655   +  656      *     ° 


6       +6       +6      +6       -:-  t         »     0 
60  63  64  65  56 
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i.e. 


(6°    h   (cos  a   +T-A^'  +  (Fn  "  S£T)xi    +  (Fi2  "  ^4^2 
10  d      Ebd  11       Ebd      1  12       „^      ■: 


2  2 

+  ?£_x     +i£_X/     =    o 

Ebd    3        Ebc,2    4 


2  2  2 

(4  "  ~2-K\  >   +   (F9,    "  ffT)Xl     +   (F22   +  — 3)X2 
20       Ebd       aD  21        tbd       '"  Z2        EbdJ       ' 


Ebd2    3         Ebd3    4 


(5.3) 


%  +   (!'ina    +T-)-EbT-)   +EbdXl     "  ^7X2    +   (F33   "  ESI>X3 


+  (F34+fj)X4    ^   F35X5+F36X6     =     ° 

Ebd 


,  w         l2r  ew    .     ,        6r2  12r  ,  6r  „     .„ 

(6,      + ;r-K       )   +  5-X.     -  rX       +   ( r-  +  F      )X 

40  m3      aw  rM2    1         ,,.  .3    2  2  43       3 


Ebd 


Ebd  Ebd  Ebd 


+  <FW+r3>X4    +  F45X5   +  F46X6 

Ebd 


6,-  +  F_,X,   +  F     X     +  FX     +  F     X        =     0 

50  53   3  54  4  5j   5  56   6 


6       +FX+FX+FX+FX        =0 

60  63   3  64   4  65   5  66   6 
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Fig.  5-3  Dimensions  and  the  relation  of  the  dome-ring  and  the 
conical  shell  wall 


According  the  dimension  shovn  in  the  Fig.  5-3,  and  the  data  already 
obtained  from  previous  calculations,  those  equations  can  be  solved.   The 
coefficient  matrix  A,  A   and  the  constant  term  matrix  see  table  5.1,  5.2 
and  5.3.   The  X  value  is  listed  in  table  5.4. 


TABLE    5.1    CCtF.     MATRIX    A 

22630.19000000  1149.0861    )i  ' 

1421.              I  O.OC 

46318.624    3i  933.4162 

60(  .         ■  0000  O.uOOi  0000 

163'.  .00001  000  -1421  .00.00<  . 

1421.00000000  .0^173730 

1421.00000000  i-ffOC  •     - 

799.98.'  -.00015800 

.  I  0  C  JOu  O.UOOOOOOO 

30.37956500  -.00173757 

■'.  000  o.ooc  ooooo 

"10.51054700  .J4182137 

TABLE  5.?  INVERSION  OF  MATRIX  A 

-.■  <■'  (  2348  .00003257 

.   '  I  3779  .001  1286' 

.       i  58474  -. 00026229 

-.       i  B4751  -.U0158927 

-.    001  247b  .00000786 

•UU043882  .Oul6il66 

.00063327  -.U0031556 

-.-,.'.41992  -.wu265737 

•u386550^  -.01926180 

-.02563610  -1.161121  . 

-.00297175  .ool4815u 

■  '  190974  .16717563 
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1630.00000000 

0.0000 

-1421.00001 

0.00000000 

-1629.07200000 

23.93204600 

1420.48800000 

-10.017272CO 

.92809668 

6.44339980 

.32085286 

2.20394340 

-.00001887 

.00000077 

.00007396 

-.00000885 

-.00026192 

.00000993 

.00054121 

-.00001645 

.03303850 

2.91942280 

-.00250509 

.°0086341 
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TABLE  5.3  CONST   MATRIX 
15(  4189.  ■ 

-16102173. OC 

-5&&57982.O01  00000 

7325131.50000000 

-38468895.'  000000 
124701. 5900C  I 

ft  <:LE  5.4  X  FORCES 


Xt     55840. 697U000C 
X.    -45645.55300000 
Xl     74517.989-' 
X(  43339679.00000000 
Xt  -6303357. OuOC 
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CONCLUSION  AND  DISCUSS  IOHS 

After  the  edge  forces  X  are  obtained,  the  corresponding  forces  produced 

in  the  shell  due  to  those  X's  can  be  achieved.  Let  &   represents  the 

i 

corresponding  force,  N  ,  N  ,  M.,  M  ,  etc.,  of  the  shell,  H   represents  the 
0   8   0   0  i 

force  due  to  unit  corresponding  displacement  which  was  shown  in  table  (3.7). 
Thus 

Cy.      =  F.X.H.  (5.4) 

Combined  the  (i.  and  the  relative  force  from  the  membrane  theory,  the 
resultant  force  is  obtained.   In  Fig.  6-1  are  shown  the  forces  at  edge  of 
the  shell  both  due  to  membrane  theory  and  due  to  the  edge  force  X.  Adding 
those  two  forces  together,  the  total  resultant  forces  can  be  obtained. 

From  the  membrane  theory  the  ring  is  in  tension  and  the  dome  is 
entirely  in  compression;  and  the  hoop  force  in  the  shell  wall  is  in  tension 
while  the  meridian  force  is  in  compression.  The  membrane  ring  tension  is 

T  .  HQ  =  N^rCosa 

i.e. 

T  r-_   (11331.6  +  15073.5)  Cos  30°  x  100  =  2,286,679.66  lb 

This  force  is  resisted  by  the  steel  alone.  Allowing  a  tensile  stress  of 
20,000  psi  In  steel,  then  114.3  sq.in  of  steel  is  required.   That  means 
29  -  ?1£  bars  have  to  be  provided.   However  in  designing  a  large  ring 
protressing  should  always  to  be  used.  Assume  the  final  stress  of  pretressed 
steel  is  of  120,000  psi 
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I 


i 


u 
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2,286,679.65     ,„  ,     , 
120,000   "  "  L  sq-ln 


Which  may  be  supplied  by  25  wires,  each  1  inch  in  diameter.  An  initial 
stress  of  150,000  psi  required  to  compensate  for  the  assumed  losses  of  20 
per  cent.   In  Fig.  6-1  also  show  that  the  action  of  the  ring  in  collection 
the  tension  forces.  An  increase  in  ring  size  and  corresponding  decrease 
the  dome  and  wall  hoop  tension  and  bending  moments. 

The  reinforcement  arrangement  for  both  the  dome  and  the  conical  shell 
wall  not  only  takes  care  the  direct  tension  or  compression  stress,  but  also 
must  provided  the  among  of  steel  to  resist  the  bending  moment  near  the 
junction  of  it.   Hence  it  is  beyond  the  scope  of  this  report,  there  is  no 
further  discussion  here. 

Another  important  factor  for  stress  analysis  in  the  shell  is  the  effect 

of  the  volume  changes  due  to  temperature.   Suppose  the  temperature  variation 

o 
is  100  F,  the  radial  displacement 


H  =  .  fe  tra  =  -  0.0000063  x  100  x  100  a  0.065  ft 


The  correspondent  moment,  meridian  and  hoop  force  due  to  this  temperature 
change  can  be  obtained  simply  from  multplying  this  displacement  with  the 
value,  by  means  of  stiffness,  which  have  already  been  obtained  in  previous 
calculations* 

It  is  sometimes  important  to  include  the  possible  effects  of  wind  load, 
eqrthqur.kes  and  blast  loads.   In  analysis  of  the  forces  due  to  such  loadings, 
generally,  a  loading  pattern  is  assumed  which  can  be  expressed  mathematically. 
All  of  these  observations  should  be  considered  in  desiging  of  shells.   It  is 
not  intended  that  specific  values  given  here  be  applied  on  the  structure  in 
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this  report. 

In  analysis  of  the  axisymraetric  shells  by  the  use  of  membrane  theory, 
it  is  found  to  be  quite  simple.   But  generally  this  method  is  exact  only 
for  shells  In  which  the  practicular  solution  to  the  shell  equations  coin- 
cides with  the  membrane  solution.   In  other  words  if  the  surface  load 
produces  bend  1713  of  the  shell  the  membrane  theory  is  not  accurate  and  may 
lead  to  great  error  in  the  solution.   The  analysis  presented  in  this  report 
combines  both  the  membrane  theory  and  bending  analysis.   Thought  it  is  an 
approximate  method  in  using  finite  difference  method  to  solve  the  differ- 
ential equations,  yet  till  now  it  is  the  only  method  that  can  be  applied 
in  solving  the  parabolical  shells.   The  numerical  calculation  in  this  report 
merely  illustrate  the  principle  involved.   It  is  found  that  the  solution  is 
reliable  and  satisfactory. 
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NOTATIONS  AND  ABBREVIATIONS 
Notation. 

Et3 
n       _     _     Bending  rigidity  of   shell. 

12    (1   -V") 

E   a  Young' fi  modulus. 

F   =  Plaxbility. 

H   =  Stiffness. 

k   -  Constant  defining  the  shape  of  paraboloid. 

M,  -  Meridional  bending  moment  due  to  membrane  theory. 

M   -  Tangential  bending  moment  due  to  membrane  theory. 
8 

Mj  ~  Meridional  bending  moment. 

M  =  Tangential  bending. 

Nj  =  Meridional  membrane  force. 

N  =  Hoop  membrane  force, 

e 

V.     =  External  load  per  unit  area  acting  in  the  direction  of  the  tangent  to 

the  meridian. 
P   -  External  load  per  unit  area  acting  inwards  normal  to  the  shell. 
p„   -  External  load  per  unit  area  in  the  direction  tantent  to  the  circular 

sect  ioK. 
p   =  Intensity  of  horizontal  projecting  load. 
q   =   Intensity  of  uniform  distributed  load. 
Q   k  Shear  force, 

r   -  Radial  distance  from  axis  to  point  on  the  shell. 
r,   =  Meridional  radius  of  curvature. 


98 


r,  -  Circumferential  radius  of  curvature. 

R  =  End  load. 

s  =  Distance  along  the  shell  in  the  meridional  direction. 

t  =  Thickness  of  the  shell. 

U  =  r2Q 


ian. 


v  =  Angle  of  rotation  of  a  tangent  to  a  meridis 

w  =  Displacement  normal  to  shell. 

l)  =  Poisson's  ratio. 

AH  =  Radial  displacement. 

60  =  Circumferential  angle  of  meridian. 

0  =  Angle  between  axis  of  shell  and  normal  to  shell. 

K  =  Modulus  of  fundation. 

60  =  Strain  in  meridional  direction. 

g  -     Strain  in  horizontal  direction. 
0 

Abbreviations. 

psi   =  Pounds  per  sq.  in. 

lb   =  Pound (s) 

in   a  Inches. 

ft   =  Feet  , 

sq   =  Squre 
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C   C   PROPERTY  CF  THE  DOME 

DIMENSION  A(11»6).B(11.5I 

PI=3. 1415926538 

DA=PI/60. 

A<1.1)=0.(, 

DC  1C  1=2,11 
J  =  I-1 
10  «(Iil|  =  MJ,l  )+na 
DC  20  1=1,11 
A<I.2)=SINF(A(I,1)> 
A( I ,3)=C0SF( A( 1,1)) 
A( I ,4)=1./A( I ,3) 
2   A( I,5)=A< I,2)/A( 1,3) 
All ,6)=  . 
DA=3.C 

DO  30  1=2,11 
J  =  I  -1 
30  A ( I  ,6)=A( J.6  1+DA 
PUMCH  6 

"UNCH  b,( (A( I,J),J=l,6) ,1=1,11) 
CK=SORTF(3.)/60O. 
DC  40  1=1,11 
6 ( I  ,1 ) =A( I  ,6 ) 
P(  I  ,3)=A(  I  ,4)/(2.*C<) 
6(  I  ,?)=B(  I ,3 )*A( I ,4) K#2 
H( 1  ,4)=R(  I,3)*A( 1 ,2) 
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40  B( I  i5)=CK*B( 1 >4)**2 
PUNCH  7 
PUNCH  0»( (B< I .J) tJ=l »S )  »I=1 »11 I 

5  -CRMAT(bF12.6,F6.2) 

6  FCRMAT(6X,3HRAD.7X>4HSINX»8X,4HCCSX,8X,4HSECX,8X,4HTANXt6Xi3HDEG/> 

7  FORMAT </MX.3HDEG»7X.2HRl»13X.2HR2tl4XilHX.14X.lHY/ I 

8  FCRMAT(F8.2.4F]4.b ) 
STOP 

END 
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C  MEMBRANE  FORCE  (DEAD  LOAD) 

DIMENSION  Alll.6liBlll.5!.CI11.5l 

RFAD,((A(I,J),J=1,6),I=1»11) 

READ.I  <B(  I ,J) »J«1»5)  .1  =  1.11) 

CK=SQRT(3. 1/600. 

C  =  75. 

PI=3. 14159265 

1=1 

C  (  I  .1  )=■  . 

C(  I  .21=  >.l 

C( I »3)=-Q/(4.»CK) 

C(  I  »4)=C<  1.3) 

DC  10  1=2.11 

J=I 

C(  I  .1 )=b( J.l ) 

W=1.+4.*CK**2*B< I .4) **2 

S=SQRT(W) 

C( I . ?)=PI*0*(S**3-1. )/ (6.*CK**2 ) 

C( I .3)=-C( I,2)*S/(4.*PI*CK*R( J,4)**2) 
18  C(I,4)=C<I .21/(2. *P!*B( J,2)*A ( J.2)**2) 
10  C(  I  ,4)=C(  I  ,4)-CS*fi(  J»4)*A(J»3)/A(J,2) 

3UNCH  15 
17  PUNCH  16. (  (C(  I ,J) .J  =  l  .4) ,1=1,11  ) 

15  FORMAT  (//4X,3HDEG,13X,1HR,17X,2HNM,16X,2HNR) 

16  FORMAT  (F8.2.3F18.4) 
STOP 

END 
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C   C   MEMBRANE  FORCES  DUE  TC  SNOW  LOAD 
DIMENSION  A(  11  .6)  ,P( 11,5) .01 11.4) 
RFAD, (  (A(  I  ,J)  ,j= 1,6)  ,1  =  1,11 ) 
READ, I (B( I, J) »J»1»51 ,1=1,11) 
CK=S0RT(3. ) /6OO4 
PJ=30. 

PI=3. 14159265 
1  =  1 

D  <  1 , 1 1  =  0  .  (. 
D(  I  ,2)=  ... 
D(  I  f3)=-PJ/(4.*CK) 
D( I  ,4)=D(  I ,3) 
DC  K  1=2,11 
D( I ,1 )=R(  I  ,]  ) 
W=1.+4.*CK**2*B( I ,4  1**2 
S=SORT(W) 

D(  I  ,2)=PI*PJ*l-  (1,41**2 
D(I,3)=-D(1,2)/(2.*P1*B(I,41*A(I,2  1) 
1    MI,4)=D(I,2I/(2.*PI*B{I,2)*A(I,2)**2I-PJ*d(I,3)*A(I,3)**2 
PUNCH  15 
°UNCH  16, ( (C(  1  ,J) , J=l,4) ,1  =  1 ,11 ) 

15  FORMAT  ( //4X.3HDEu,13X,lriR,17X,3MNf-iP,16X,3HNRP) 

16  FORMAT  (F8.2,3Flb.4) 
STOP 

END 
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IC   DISPLACEMENT  FROM  <-     THEORY 

D I ME NS I C N  A(ll,6)tB(lli5)tC<lI.4),D(ll,4),G(llil3) 

READ,( (A( I,  J)  ,J=1, 61.1  =  1,11) 

READ»((R(!,J),J=1,5),I=1,11) 

REAU,( (C(  I .J) »J«1»41  ,1  =  1,11  ) 

V  «  0  .  2 

T»0.5 

Gll.ll  )  =  .  ., 
GUtl2)«0.( 
0(1  ,131=0,1 
DC  If.  1=2,11 
9  G( I ,1)=B(  I  ,1 ) 

5( I  ,2)=B(  I ,3)*A(  1 ,2)/T 
G(I  ,31  =  1 ,/(b(  I ,2)*T*A<  1,5) I 

G(  1  ,A)=G(  I  ,3)*(C(  1  ,3  )*([(<  I  ,2)+v*h(  i  ,3)  )-C(  I  ,i,l*(B<  I  ,3)+V*b(  1,2))) 
0( I ,5)=-2./(A< I ,2)**2*A( I ,2) ) +A ( I ,4)**2/A( I ,2)+2.*A( I ,3)/A( I ,c  | 
0( 1  ,6)=2./(A(  I .5)**2*A< I .5  )  )-3.*A(  1 ,4>*a( 1,5) 
G ( I  ,7)=G(  I ,b )+G( I ,6) 

G( I ,6)=-A( I ,4)**2/< A( I ,2)*A( I ,5) 1+3. *A< 1 ,4)**2*A< I ,4)**2 
G<I,9)=A(I,4)/(A< I,2)#*2*A( I ,2) ) -A ( 1 ,4) /A (I ,2) 
G( I  ,10)=0( I , 7)-v*(G(  I,8)+G(I,9)+1./A(1 ,2)**2*A( 1,3)) 
G(  I  ,11)=B( I ,1 ) 

G(I,12)=G(I,2)*(C(I,4)-V*C(I,3)) 
10  G(  I  ,13)  =  (G(  I  ,4)-G(  I,lc.  )  )/b(  I  ,2) 
P!">'CH  5 
PUNCH  6,((G(I,J),J=ll,13),l=i,ii) 


105 


5  FORMAT  (//4X.3HDEG.8X»6HE*D(H) ,1*X»7HE*U(AN) 

6  FORMAT  (F8.2.2F18.6) 
STCP 

END 
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C   C   CCEFFECIENT  MATRIX  CF  AX=G 
DIMENSION  A( 11.6) .BI22.22) 
READ.) IA( I, J] , J=l,6) ,1=1,11) 
DA=3. 1415926536/60. 
DD=DA**2 
DC  10  1=1,22 
DC  10  J=1.22 
D4=4.*DD 

CK=S0RT(3. )/r   . 
G=6.*<l.-U. 25**2  )*DD/CK 
10  B( I,J)=0. 

b(l.l)=-(2.+D0/(A(2,2)*A(2,3))**2) 

B( 1.2)=-0.5*<DD*AI 2,4)**5)/(2.*CK) 

e<1.3)=l.+DA*(-2.*AI2.5)+l./A<2.5> 1/2. 

B(2,l)=G*A(2,4)**5/( ^.5**3 ) 

b(2.2)=b< 1,1 ) 

b(2,4) =  B( 1 ,3 ) 

DC  20  1=3,19,2 

J  =  l+1 

K=I  +  2 

L=  I  +3 

N*( 1+3)/? 

M=I-1 

KL=I-2 

IMI,KL)=1.+(2.*A(N,5)-1./A(N,5) )*DA/2. 

b(  I  ,1 )=-(2.  +  DD/(A(N,2)*A(N,3)  ) **2 ) 
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o(  I  ,J)=-u.5*DD*A(N»4 1**5/1 2.*CH) 
B(I tK)=l.+DA*<-2.*A(N»5)+l./A<N»5) 1/2. 
B(J»M)=B( I »<L) 
B(J»I)=(G*A(N» A) **5 1/0.125 
B(J»J)=B( 1 •! I 
20    B<J»L)=BI I ,K> 

DK=A( 1 1,3 )**2*A( 11,21/ (2.*0.5*DA) 
DC    30     1=1,22 
DO    3u    J=l,22 
BI21»17)=DK 

b(21il9)  =  (J.2/(0.5*A(  11,2)1 
b(21,21 >  =  -DK 
i  (22,2^  )  =  l.i. 

PUNCH  6,( (t( I ,J) ,J=1 ,22) ,1=1,22 ) 
6  FORMAT  (4F16.fi) 
STOP 
END 
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C   C   INVERSION  CF  MATRIX  A 

DIMENSION  AI22.22) »B(22) >C!22) 
12  FORI-. AT  13E2^.S) 
READ.N 
NN=N-1 
11  FORMAT  I3F2A.0) 
10  READ»( (A( IfJ)»J*l»N) »1«1»N) 

PUNCH  5»{(A(I.J)»J=1»15).I=1»15) 
5  FORMAT  (  1GF7.1 ) 
A(l»l)=1.0/A(l»l) 
DO  110  M=1,NN 
t«M+l 
50  DO  6u  1=1 »M 
B(I)=0. 
DO  60  J=1,M 
60  B( I )=p ( I )+A( I ,J)*M J,KI 
D  =  0.0 

DO  70  1=) iM 
70  D=D+A(K, I )*E ! I ! 
D=-D+A (KiK) 
A(K»K)=1./D 
DO  80  I=1,M 
8.'  A(  I  iK)=-B(  I  )*A(K>K) 
DO  9u  J=ltM 
C(J)=0. 
DO  90  1  =  1  tM 
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9C  C(J)=C(  J  )+,'(<.!  )*A(  I  ,J> 
PC  100  J=1,M 

]    A<K,J)=-C(J)*A(K»K] 

DO  110  1=1, M 

DC  110  J=1,M 

11   A( I  ,J)=A(  I  ,J)-b(  i )*A(K.J) 

DC  170  1=1, N 
170  PUNCH  12»(A( I,J)»J=1,N) 

STOP 

END 
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C   C   FORCES  nUE  TO  BOUNDARY  DISP.  < D(H ) =1 ,D ( AN ) =0  AND  D ( H ) =0 ,D ( AN ) =1 ) 

DIMENSION  A<11.6)»B(11»5)>C(12).D(11,3> 

RFAD, ( (A( I ,J)  ,J  =  1  .61,1  =  1,11) 

READ»( (B(I.J) »J«lt5> tl«l.ll) 

C(  1  )=  . 

DO  16  LM=1,2 
11  READ, (C(  I  ) , 1=2,12  ) 

DA=3. 1415926538/60. 

D(l »1)=0.C 

D(1.2)=0.( 

D<1.3)=0.0 

DO  1C  1=2,11 

M=I-1 

L=I  +  1 

D I  I  , 1 ) =R (  I  ,  1  ) 

D( I  ,?)=-C(  I  ) /(B(  I  ,3)*A(  I  ,5  )  ) 
10  D( I ,3 )=-(C IL)-C(M) ) / I2.*DA*S (1,2)) 

PUNCH  b 

PUNCH  6,( (D( I ,J) ,J=1,3) ,1=1,11) 
.6  CONTINUE 

5  FORMAT  [//4X»3HDE6»15X»4HN(M) »15X»5HN(AN) 

6  FORMAT  (F3.2.2F20.8) 
STOP 

END 
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C   C   MOMENT  DUE  TC  BOUNDARY  DISP.  (D<H ) =1. ,0 ( AN) =0 ,AND  D (H) =0,D( AN ) =1 . ) 
DIMENSION  A(  11.6)  >B(  11.5)  »C(  12)  »Ct  11  ,5) 
READ,(IA(I.J).J=1.6).I=1.11) 
READ.I  (B(  I  ,J) »J  =  3 .51.1=1.11) 
C(1)=0.C 
DO  16  LM=1.2 
] l  READ12»(C( I ) .1=?.] 2) 
DA*3. 1415926536/60. 
V»i  .? 
D0=-0.5**3./(12.*ti«-0.2*»2) ) 

DO  1C  1=2.11 

L  =  I  +1 

M=  I  - 1 

b(  1  .31=  • 

D(1»4)=0.U 

D(l»5>=0«( 

D(  1  ,1)=(C(L)-C(M)  )/(2.*DA*B( I .2) ) 

D(I,2)=V*C( I  )/(A( 1  ,5)*FU  1,3)1 

D<I»3)=Bt!»ai 

D(  I  ,A)=DD*(D(  1  ,2)+IM  I  .1  )  ) 
10  D( I »5)=D0*(D( 1 ,2)/V+V*D( 1,1)) 

PUNCH  5 

PUNCH  6»(  (D<  1  .J) »J  =  3,5) ,1  =  1,11) 

5  FORMAT  (4X.3HDEG.13X»4HM(M>  ,  17X  ,  5mm  I  AN ) 

6  FORMAT  IF8.2.2F20.8) 
16  CONTI 
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12    FORMAT    (E24.b) 
STOP 
END 
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C   l|        FORChS  DUE  TC  DEAD  LO 

DIMENSION  A  (11.3) 

PI=3. 1415926538 

DA  =  5. 

DP=115.45 

DL=65.45 

DQ-10O. 

AA=CCSF( 1.04719744) 

AC=SINF ( 1.. 4719744) 

AB=SINF ( 1 .04719744 )/S INF (1.04719744) 

UC    10  1=2.11 

J-I-l 

All.  11=65.45 
1   A ( I ,] )=A(J,1 )+DA 

DC  20  1=1,11 

A ( I ,2)=-DQ*lDP**2-A( I .1 1**2 )/(2.*A( I ,1 )*AA) 
2^  ACI»3)=t>Q*A(  I  tl)*AB*AC 

PUNCH  5 

PUNCH  6.1  (A(  I .J) »J  =  1»3> ,1  =  1  ,1  1 ) 

5  FORMAT  ( /4X.lriY»ldX.4HN(Y) ,18X,5riN ( AN) 1 

6  FORMAT  (F8.2.2F24.8) 
STCP 

I 
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C   C   MEMBRANE  FORCES  DUE  TO  DOME  LOAD  P 
DIMENSION  A(ll»3) 
PI=3*1415926538 
Rs-2.*PI*<10Q.+l«75f#3.5*3.5*150.-3559931. 

DA=5. 

A  t 1  , 1 1=65.45 

AA=CCSF( 1.04719744) 

AC=SINF< 1.^4719744) 

=  SINF  I  1  .  ■•  4  719744  1/SINFl  1.04719744) 

DC  10  1=2,11 

J=  1-1 
10  A(  I  ,1)=A( J,1)+DA 

DC  20  1=1,11 

A(I,2)=R/(PI*A(1,1 ) #2.* AC* A A ) 
20  All »3>=-R*AB/(2.*Pl ) 

PUNCH  5 

PUNCH  6, (  (A(  I  ,J) ,J  =  ]  ,3) ,1  =  1 ,]  1 ) 

5  FORMAT  ( /4X,1HY,18X,4HN( Y) ,18X,5HN (AN) ) 

6  FORMAT  (F8.2.2F24.8) 
STOP 

END 
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C   C   DISPLACEMENT  AND  ROTATION  DUt  TC  MEMEBRANE  THEORY 
DIMENSION  A(  14,4)  ,61  Hi*) 
r=C .6666667 
G  =  1  ■-  '  . 

PI=3. 1415926538 
DA  =  5. 
D|_  =  65.45 
V  =  u.2 

R«+2.*PI» «XOO. +1.75 1*3. 5*3. 5»15G. +3559931. 
AA=COSF( 1.04719744) 
AB=bINF( 1.04719744) 
AC  =  AA/At» 

DO  10  1=2,11 
J  =  i  -1 

All, 11=65.45 
10  A (  I  , 1)=A( J,1)+DA 

DO  20  1=1,11 

AD=0/T 

/\F  =  AB*AC*A(  1,1) 

AF=1./(2.*A( I ,1 )*AA) 

AtI»2)=AD*A< I,i1*A6*(AE+AF*V*( 1 15. 5**2 -A ( 1,11**2 ) ) 
?_•      All  ,3)=AD*AC*AF*(A(  I  ,1)*(A(  1  ,1  )  -V  )  -1  1  5  .  45**2*  (  l.+V)  ) 

PUNCH  6 

PUNCH  7,1  (A(  I  ,J> ,J  =  1 ,3  1  ,1  =  1,11) 

DO  30  1=1,11 

b(  I  ,  11=A(  I ,1 ) 
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B  (  I  ,2  )  =R*  A  (  I  .  1  >*Ab*  I  AC+V/  (Atl.11  *Ad*AA  )  )  /  (  2.*P1*T  ) 
B ( I ,3 ) =R*AC* ( ( V-l • ) / ( A ( I .1 i *AA*AB ) -V*AC ) 

3   B(  1  ,3)=B(  I  ,31/(2. *PI*T) 

PUNCH  6 

PUNCH  7,(  (B(  I, J) ,J=1»3)  ,1  =  1,11) 

6  FCp:«T  (/AX,1HY,16X,4HF*DH,18X,4HF*DA) 

7  FORMAT  (F8.2,2F24.8) 
STOP 

E  ID 
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C   C   COEFFICIENT  MATRIX  OF  AX=G 

DIMENSION  M 15»15)*B(151 

V«0.2 

PI=3. 1415926538 

COT=COSF( 1.04719744) /S I NF< 1.04719744) 

uA  =  5. 

DB*0. 66666  1 

b( 1 1=55.45 

DO  10  1=2,15 

J=  I  -1 
10  B(I)=BU)+C 

DO  20  1=1,15 

DO  20  j=l,15 
20  All  ,J)='  .1 

A(2,2)=-l. 

A(2,4)=i. 

A  (  3  ,  3  )  =  1  .  0 

A(14,12)=-l.u 

A (14, 141=1. u 

A(15»13)=l.t 

DO  30  1=3,13 

A ( 1 ,1 ) =(B( 3) /DA-1. 1 /DA**3 

A  !  1  ,2  1  =  (-4.*BI  3  )  /OA  +  2  .  1  /DA**3 

A(],?)=(9(3)/nA**4+2.»<CCT)**2*( l.~V**2 1 /(DB**2*B( 3 1 1 1 *6. 

A (1 ,4)=(2.*P(3)/DA+1. )*(-2. ) /DA«*3 
31  A  (  1,5  )  =  (  B(3)  /DA+1.  1  /DA**3 
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DC   40    1=4.13 
N=  1  +  1 
R=l+2 
M  =  I  - 1 
L=I-2 

A(I.L)=(B<I) /DA-1 . )/DA*«3 
ft(I,M)»(-i.*B<I»/DA+2.)/OA«»3 

A(I,I)=(B(I) /DA**4+2.*(CCT)**2*(1«-V**Z )/(DB**2»b( I ) ) 1*6. 
A(  I,N)=(2.*BU)/DA+l.)*(-2.>/DA**3 
40    A<  I  ,<)  =  <!<(  I  )  /DA  +  1.)/DA**3 

PUNCH  5t( (A(l*V)»J«lfl5)»I"lfl5) 
5  FCRMATI4F18.6] 
STCP 
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C   C   CALCULATION  CF  Z-VECTCR 

DIMENSION  A(15,15),B(22),C(20»5) 

12  FORMAT  (3E24.0) 

DC  15  1=1,15 
15  READ  12, (All  ,J)  ,  J  =  1  , 1  5  ) 

C(] ,1 1=55.45 

DA  =  5. 

DO  180  1=2,15 

J*l-1 
18  J  C( I  ,1)=C( J,l )+DA 

DO  200  K*l»4 

READ,(d(  11,1  =  1.15) 

DO  120  1=1,15 

Cf I »2)=C.CJ 

DO  120  J=l,15 
120  C( I ,2)=C( I »2)+A( I ,J)*B( J) 

PUNCH  130 
13   FORMAT  (/4X,lHY,12X,aHZ-VECTCR> 
14:  PUNCH  6,( (C( I, J) ,J=1,2) ,1=1,151 

6  FORMAT  (F8.2.F24.B) 
200  CONTINUE 
ST~P 
END 
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C   C   FORCES  DUE  TO  EDGE  EFFECT 

DIMENSION  A( 15.2 >  »R( 15.4) »C< 15.3) 

TAN  =  SINF(  1.04719744)  /CCSFI  1.04719  7  44) 

V  =  l  .2 

T  =  ^ .666067 

DA  =  5. 

K=T**3/ (i 2. *(1.-V**2)*CCSF< 1.04 7 1944  1**3) 

DO  200  K=l»4 

READ»( ( A( I >J) > J=1.2 ) » 1=1 » 151 

DO  10  1=3.13 

L=I-2 

M=I-1 

M  =  I  +  1 

LM=I+2 

AA=-A(L»2)*A(  I,1)/DA**3  +  A(N'.2)*(A(  1*1  )  /DA+1  .  )  /DA**  2 

AA=AA-A( I,2)/PA**2+A(N.2)*< l.-AU , 1 ) /DA ) /DA**2 

AA=AA+A(LM.2)*A< I . 1 ) / ( 2 . *DA**3 ) 

l  ( I ,1)=A( I  ,1) 

Bl I ,2)=k*AA 

d( I .3)=b( I >2)*TAN 
10  B(I»4>=T*A(I»2)/(A(r.l)*TAN! 

PUNCH  5 

PUNCH  7il  (B(  1  .J) .J=1.4) ,1  =  3.13) 

DO  20  1=3.13 

C (  I  ,1 )=A(  I  ,1 ) 

C  (  1.2  )=-!?*(  AIM.  2  )-2.*A(  I,2)+A(N,2)  ) /DA**2 
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2   C( I ,3)=V*C< I  ,2) 
PUNCH  4 

PUNCH  6> ( (C( I »J) »J=1 .3 ) ,1=3.13) 
h    FORMAT  (//4X»1HY»18X»4HM(Y) »20X.5HM(AN) I 

5  FORMAT  ( //4X,lHY»16X»lHCw,l&X,4hi\|(  Y  )  >  !6X,bHN(  AN) 

6  FORMAT  (F8.2.2F24.8) 

7  FORMAT  (F8.2.3F2U.8) 
CONTINUE 

STOP 
END 
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C   C   CALCULATION'  CF  FLIXBILITY  MATRIX  (  n'ALL  ) 
DIMENSION  AI22.22) »B<22) tC<22).E<22) 

11  FORMAT     13F24.8) 

12  FORMAT     (3E24.6) 
RFAD.N 

NN=N-1 
DO  10    1=1.4 
AA=SINF<  1.1.4719744) 
A(1»1)=-0.53652068*AA 
A(1,2)=1.2812699*AA 
Atl.3)=0.0C 194688*AA 
All »4)*0.0(  332947*AA 
A  (2.1)='  .   --15366 
A<2>2)=-0.  430242 
A(2»3)=0. 03292389 
A (2. 4  I =  .    I  094 
A (3 tl)=-0. 00165074* AA 
A(3»2)=-0.0o460'3333*AA 
A(3,3)  =  1.17l7cS72*AA 
A (3t4)=3.3840397*AA 
A(4,l  )=-',  ...2376767 
A  <4»2)=-0'.  0430242 
A (4,3 ) =  .      257 
1<  A(4»4)=0.(   10094 
All  ,1  )  =  ]..<  /A  (,1.1) 
DO  110  M=1,NN 
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<  =  !■.+  ! 

5C    DC    6C     [=>1»M 
B  ( I )  =1  • 

DC    60    J=] >M 
60    B(l)=BI  1>+A(I.»J)*ACJ»K) 

D  =  '  . 

DC    TO     1  =  1  ,1' 
7      D»D+A< K»l)*8t I) 
D=-D+A(K.KI 
A(K.K)=1./D 
DC    80    I  =  l.i- 
81    A(I.K)=-B(1)*A(K»KI 
DC    ?Q    J-l  »  ■■ 
C(J)=     . 
DC    90     1  =  1  iC 
90    C( J)=C( Jl+AIK.I 1*A(  I  .J) 
DC    100    J=1.M 
10  0    A(K.J)=-C( J)*A(K.< 1 
DC    110     I  =  1.M 
DC    110    J=1»M 
110    A( I tJ)=A( I >JI-H( I >*AU.J) 

DC    170     I»1»N 
17-     PUNCH    11,(A(1»J)»J=1»N) 
DC    2o    1=1.2 

£  1 1 )=15^.**l.b*33.*SQRTF(40uC. ) 
20    E  (2  1=151  .**  1.6*33.  *SURTF(40t.C.  1*144. 
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PUNCH  5 

PUNCH  6»(E(I)iI«l»2) 

5  FORMAT  (/MX,7HE-VALU£) 

6  FORMAT  (2F24.6) 
STOP 

END 
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C   C   CALCULATION  OF  X  VALUE 

DIMENSION  A(22,22)  ,B(22),C(20*2)>tUl) 

11  FORMAT  (3F24.6I 

12  FORMAT  (3E24.6) 
READ.N 

I-  v  =  N-l 
1   READ,< (A(I»J)»J"1»N)»I=1»N) 

PUNCH  15 
15  FORMAT  t7/4X,3H(l)»lX.4HCCEF.lX.6HMATRlX.lX»lHA) 

PUNCH  1 1 , ( (A( I .J) ,J=1,N> .1=1 iN! 

All  ,1  1  =  1.-- /A  (  1,1) 

DO  U(  [•'=], NN 

K»M+'l 
50  DO  60  I  =  l,i' 

B  ( I ) =0 . 

DO  6w  J=1.M 
6  J  B  i  I  )=D(  1  )+A(  I  ,J)*A(  J.K.) 

i  ■  =  . 

DO  70  1=1. M 
70  D=D+A(K, 1 )*B( I ) 

D=-D+A(K,K) 

A (K,K)=1./D 

DO  80  !=!,■ 
8u  All  «K)=-B<  I  )*A(K.K> 

DO  90  J=1,M 

C ( J  )  =0. 
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DC  9C  1=1 fM 
90  C(J.)=C(J>+A(KtI)*A<  I»J) 

DC  100  J=1,M 
IOC  AK.J)=-C(  J)*A(K.K) 
DC  110  I  =  1,M 
DC  110  J=1»M 
110  A(  I  ,J)=A(  I  ,J)-B(  I  »*A(K.JI 
PUNCH  14 
14  FORMAT  <  //4X,3H(2)  .1X.9HA-INVERSE! 
DC  17U  1=1. N 
17.  PUNCH  11,  (  A(  I  ,J)  .J=1.N) 
WTH=67431871.4 
WTA=-89831.86 
VBH=38468895.2 
WBA=-1 2470 1,59 
AA=815. 
AH=1421. 
Ai=800. 
DDH=-1484699. 
DDA=-49.7 
XD=-11331.6134 
XWT=-13073.541 
NWB  =  -54(,29.318 
DC  20  1=1,6 
Bl 1 )=-DDH-1.72*XD 
d(2 >=-DDA+AD*XD 
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B (3 )=-*TH-i..d77*AA*XWT 
B t4l=-WTA-AE*0.6*XWT 
8 (5)=-WBH 
20  B(6)=-WBA 

PUNCH  5 

PUNCH  6,<  B( I  ) ,1  =  1,6) 

5  FORMAT  (  //hX,1X,3h( 2 ) ,1X,5HCCNST> 1X.6HMATRIX) 

6  FORMAT  (F2A.8) 
DC  12C  1=1,6 

ci i >=o. : 

DC  120  J=l,6 
12 C  C( I )=C( 1 )+A( I ,J)+AI I ,J)*b( J) 

PUNCH  13 
1'   FORMAT  !//4X,3H(M ,lX,lHX,lX,bHFCRCt) 

PUNCH  6, ( C( I ) , 1=1,6) 

STCP 

END 
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The  purpose  of  this  report  was  to  introduce  the  design  of  thin  shells 
by  (!)  de-riving  the  membrane  equations  for  the  elastic  analysis  of  shells  of 
the  form  of  a  surface  of  revolution  and  loaded  symmetrically  with  respect  to 
the  axis  (2)  by  describing  the  physical  behavior  of  a  well  defined  system. 
In  order  to  illustrate  the  formula  derived,  the  dome  structure  is  presented. 

Due  to  membrane  theory  the  stresses  analysis  of  a  shell  structure  some- 
times is  not  true.   Hmoshenko  develops  the  general  homogenous  equations  for 
an  axisymmetric  shells.   To  solve  a  shell  problem  a  membrane  solution  is 
superimposed  upon  the  solution  of  the  homogeneous  equations.   The  complete 
procedure  in  analysis  of  a  shell  structure  can  be  outlined  as  follows: 

1.   Calculation  of  forces  due  to  membrane  theory. 

?.   Calculation  of  the  displacements  at  the  boundary  of  the  shell  from 
the  membrane  theory. 

3.  The  corrections  correspond  to  unit  edge  effects  derived  from  the 
solution  of.  the  homogeneous  equations  commonly  referred  to  as  the  rending 
theory. 

4,  Compatibility  is  obtained  by  determining  the  size  of  the  corrections 
required  to  remove  the  errors  in  the  membrane  theory. 

In  design  a  shell  structure  not  only  consider  the  effects  due  to  the 
Symmetrical  surface  lo3d  and  dead  load,  but  also  need  to  consider  the  effect 
due  to  wind  load,  temperature  etc.   However,  it  is  beyond  the  scope  of  this 
report,  the  author  only  mentioned  ic  in  the  part  of  conclusion  and  discussions. 


